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Insect  Community  Responses  to  Climate  and  Weather 
Across  Elevation  Gradients  in  the  Sagebrush  Steppe, 
Eastern  Oregon 

David  S.  Pilliod  and  Ashley  T.  Rohde 


Executive  Summary 

In  this  study,  the  U.S.  Geological  Survey  investigated  the  use  of  insects  as  bioindicators  of 
climate  change  in  sagebrush  steppe  shrublands  and  grasslands  in  the  Upper  Columbia  Basin.  The 
research  was  conducted  in  the  Stinkingwater  and  Pueblo  mountain  ranges  in  eastern  Oregon  on  lands 
administered  by  the  Bureau  of  Land  Management. 

We  used  a  “space-for-time”  sampling  design  that  related  insect  communities  to  climate  and 
weather  along  elevation  gradients.  We  analyzed  our  insect  dataset  at  three  levels  of  organization:  (1) 
whole-community,  (2)  feeding  guilds  (detritivores,  herbivores,  nectarivores,  parasites,  and  predators), 
and  (3)  orders  within  nectarivores  (i.e.,  pollinators).  We  captured  59,517  insects  from  176  families  and 
10  orders  at  the  Pueblo  Mountains  study  area  and  1 12,305  insects  from  185  families  and  1 1  orders  at  the 
Stinkingwater  Mountains  study  area  in  2012  and  2013.  Of  all  the  individuals  captured  at  the 
Stinkingwater  Mountains  study  area,  77,688  were  from  the  family  Cecidomyiidae  (Diptera,  gall  gnats). 

We  found  that  the  composition  of  insect  communities  was  associated  with  variability  in  long¬ 
term  (30-yr)  temperature  and  interannual  fluctuations  in  temperature.  We  found  that  captures  of  certain 
fly,  bee,  moth,  and  butterfly  pollinators  were  more  strongly  associated  with  some  climate  and  vegetation 
variables  than  others.  We  found  that  timing  of  emergence,  as  measured  by  first  detection  of  families, 
was  associated  with  elevation.  When  analyzed  by  feeding  guilds,  we  found  that  all  guilds  emerged  later 
at  high  elevations  except  for  detritivores,  which  emerged  earlier  at  high  elevations.  The  abundance  of 
most  taxa  varied  through  time,  mostly  in  response  to  temperature  and  precipitation.  Of  the  pollinators, 
bees  (particularly,  Halictidae  and  Megachilidae)  peaked  in  abundance  in  late  June  and  early  July, 
whereas  butterflies  and  moths  peaked  in  August.  Flies  peaked  in  abundance  in  July. 

Overall,  our  interpretation  of  these  patterns  is  that  insect  communities  respond  positively  and 
negatively  to  weather  and  local  vegetation  more  than  to  long-term  climate.  Given  increasing  variability 
in  weather  and  high  probability  of  extreme  weather  events,  insect  communities  in  sagebrush  steppe  also 
may  experience  considerable  fluctuations  in  composition  and  abundance,  as  well  as  phenology.  These 
findings  have  implications  for  many  ecosystem  services,  including  pollination,  decomposition,  and  food 
resources  for  predatory  birds  and  other  vertebrates. 
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Introduction 


Climate,  in  large  part,  determines  the  distribution  of  plants  and  animals.  Reeent  aeeelerated 
elimate  ehange  eaused  by  human  aetivities  appears  to  be  no  exception.  Globally,  lower  snowpacks, 
earlier  springs,  new  rainfall  patterns,  and  warmer  temperatures  have  been  measured  over  the  last  50-60 
years  (Karl  and  others,  2009).  Since  the  late  1800s,  the  planet  has  warmed  approximately  0.85°C, 
mostly  because  of  an  increase  in  minimum  annual  temperatures  (Pachauri  and  Meyer,  2014). 

Populations  of  plants  and  animals  may  respond  to  changing  climate  by  dispersing  to  new  sites  with  more 
suitable  conditions  (climate  tracking)  or  adapting  genetically  or  phenotypically  to  the  new  environmental 
conditions.  However,  the  natural  expansion  and  contraction  of  species  ranges  associated  with  climate 
deviations  may  be  influenced  by  human-associated  barriers  or  stressors,  such  as  land  use,  natural 
disturbances  (like  fire),  or  invasive  species  (Vitousek,  1994;  Vitousek  and  others,  1996;  Sala  and  others, 
2000).  If  species  are  unable  to  move  or  adapt,  then  they  are  more  likely  to  go  locally  or  even  regionally 
extinct  (Hill  and  others,  2011). 

Investigations  of  range  shifts  associated  with  recent  climate  change  often  are  focused  on  the 
leading  edge  of  distributions  with  colonization  occurring  along  high  latitudes  and  elevations,  where 
environmental  variables  are  more  favorable  under  changing  conditions.  A  study  that  measured  the 
distribution  of  37  species  of  dragonflies  in  Britain  concluded  that  the  ranges  of  southern  species  are 
shifting  northward  (Hickling  and  others,  2005).  Two  species  whose  ranges  already  extended  to  the 
northernmost  limits  of  Britain  were  found  to  be  retracting  in  distribution  from  their  southern  boundary 
(Hickling  and  others,  2005).  Within  the  plant  communities  of  southern  California’s  Santa  Rosa 
Mountains,  species  were  reported  to  shift  an  average  of  almost  65  m  higher  in  elevation  from  1977  to 
2006  (Kelly  and  Goulden,  2008). 

Phenological  responses  to  climate  change  also  have  been  documented.  In  a  meta-analysis  of  61 
studies  considering  694  species  or  groups  of  species.  Root  and  others  (2003)  found  that  springtime 
events  of  some  populations  are  as  many  as  24  days  earlier  than  they  were  50-100  years  ago.  When 
several  groups  of  organisms  were  considered  separately,  invertebrates,  amphibians,  non-tree  plants,  and 
birds  had  springtime  phenologies  that  were  five  days  earlier  currently  than  historically.  Trees  had 
altered  their  phenologies  by  three  days  (Root  and  others,  2003).  It  is  unclear  whether  these  changes  are 
due  to  genetic  adaptation  or  genetic  plasticity,  but  changes  in  genotype  frequencies  have  been  linked 
with  climate  change  (Gienapp  and  others,  2007). 

Insects  are  good  model  organisms  for  measuring  the  effects  of  climate  change  because  of  their 
short  generation  times,  relatively  rapid  responses  to  disturbance,  and  large  population  sizes.  They  are 
among  the  easiest  and  least  expensive  animal  taxa  to  collect,  allowing  robust  sample  sizes  for  statistical 
analysis  using  relatively  high-level  taxonomic  identification  for  economical  and  swift  results  (Williams 
and  Gaston,  1994;  Cagnolo  and  others,  2002;  Riggins  and  others,  2009).  Moreover,  arthropods  are 
critically  important  members  of  their  communities  because  they  occupy  the  widest  variety  of  niches  and 
play  more  ecological  roles  than  any  other  group  of  animals  (Longcore,  2003).  Insects  are  among  the 
most  important  pollinators  on  the  planet  and  are  a  major  food  source  for  many  vertebrates  (McGee, 
1982;  Williams,  1984;  Johnson  and  Boyce,  1990;  Crawford  and  others,  2004). 
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The  topic  of  this  report  is  especially  timely  because  of  a  recent  Presidential  Memorandum,  the 
Pollinator  Research  Action  Plan  (The  White  House,  2015),  that  draws  national  attention  to  restoration 
and  preservation  of  habitat  for  pollinators  across  the  country.  This  plan  identifies  key  priority  research 
themes  including:  (1)  Understanding  pollinator  habitat  requirements  and  (2)  Understanding  habitat  loss, 
degradation,  and  fragmentation  effects  on  pollinators,  as  well  as  stressors  that  interact  with  and 
exacerbate  these  effects.  Additionally,  the  Greater  Sage-Grouse  Approved  Resource  Management  Plan 
Amendment  (BLM/OR/A/PL-1 5/05 1+1 792)  released  by  the  Bureau  of  Land  Management  (2015) 
identifies  priority  sage-grouse  habitat  management  areas  in  multiple  states  including  over  9.7  million 
acres  (3.9  million  hectares  [ha])  in  Oregon  and  over  2.8  million  acres  (1.1  million  ha)  in  Harney  County, 
Oregon,  where  the  research  study  areas  for  this  study  were  located.  Insects  are  a  critical  food  source  for 
juvenile  sage-grouse  (Johnson  and  Boyce,  1990;  Drut  and  others,  1994). 

In  this  study,  we  quantified  insect  community  composition  along  an  elevation  gradient  to  make 
predictions  about  relationships  between  insects  and  climate  in  a  semi-arid  shrubland.  Insects  have  been 
shown  to  experience  and  respond  to  changes  in  climate  along  elevation  gradients,  particularly 
differences  in  temperature  and  precipitation  (Warren  and  others,  1988;  Hodkinson,  2005). 

The  original  proposal  addressed  the  following  research  questions: 

•  What  is  the  current  composition  and  abundance  of  families  in  the  insect  community  across 
elevation  gradients  in  sagebrush  steppe? 

•  Does  composition  of  insect  communities  vary  between  study  areas  or  between  transects  within 
study  areas? 

•  How  does  interannual  variation  in  weather  and  habitat  affect  composition  and  abundance  in  the 
insect  community? 

•  How  do  changes  in  elevation  and  associated  habitat  factors  affect  the  phenology  (for  example 
emergence  and  diapause  [a  period  of  suspended  development,  often  during  the  winter  for 
temperate  species])  of  different  taxa? 

This  study  was  designed  in  two  phases:  a  short-term  study  of  the  differences  in  composition  and 
phenology  of  insects  along  elevation  gradients,  and  a  long-term  study  measuring  the  changes  in  species 
ranges  and  phenology  over  10  years  or  more  at  the  same  established  sampling  plots.  The  first  phase  of 
this  study  has  been  completed  and  the  results  are  reported  here.  We  organized  the  major  findings  into 
three  sections:  an  assessment  of  sampling  design  with  considerations  of  the  differences  in  measured 
variables  between  study  areas  and  between  transects  within  study  areas  (Section  I.  Assessment  of 
Sampling  Design),  an  assessment  of  the  insect  community  composition  and  associations  with  climate, 
weather,  and  habitat  (Section  II.  Insect  Community  Composition),  and  an  assessment  of  insect 
phenology  (emergence,  diapause,  and  seasonal  abundance)  in  relation  to  elevation,  weather,  and  habitat 
(Section  III.  Insect  Phenology).  The  detailed  questions  addressed  in  each  section  were  derived  from  the 
research  questions  in  the  original  proposal  and  are  listed  at  the  beginning  of  each  section. 
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Methods 

study  Area 

We  established  four  transects  along  an  elevation  gradient  within  the  Stinkingwater  (elevation 
1,150  to  1,600  m)  and  Pueblo  (elevation  1,300  to  1,850  m)  mountain  ranges  in  southeastern  Oregon 
(hgs.  1  and  2).  Study  areas  were  located  168  km  apart.  Transects  at  the  Stinkingwater  Mountains  study 
area  were  19  km  apart  and  transects  at  the  Pueblo  Mountains  study  area  were  9  km  apart.  These  four 
transects  serve  as  replicates.  Study  areas  were  within  the  jurisdiction  of  the  Burns  Field  Office,  Harney 
County,  Oregon,  managed  by  the  Bureau  of  Land  Management. 

Sagebrush  steppe  vegetation  spans  over  44.4  million  ha  of  land  in  the  western  United  States 
(Miller  and  Edelman,  2000),  and  more  than  80  percent  of  this  land  is  affected  by  anthropogenic 
activities  (West,  1999).  Reductions  in  biodiversity  due  to  climate  change  are  predicted  to  be  most  severe 
in  Mediterranean  climate  and  grassland  ecosystems  (Sala  and  others,  2000),  such  as  sagebrush  steppe. 
These  systems  may  be  particularly  susceptible  to  climate  change  because  of  their  aridity,  and  high 
amount  of  habitat  conversion  and  low  amount  of  protected  land  relative  to  other  major  habitat  types 
(Hoekstra  and  others,  2005). 

In  addition  to  climate  change,  land  use  and  non-native  species  are  predicted  to  effect  semi-arid 
grasslands  and  sagebrush  steppe  (Sala  and  others,  2000).  Historical  mismanagement  of  cattle  grazing, 
invasive  annual  grasses,  and  fire  regime  changes  cause  fragmentation  and  reduction  of  habitat  that 
complicate  organisms’  responses  to  climate  stressors  (Knick,  1999;  fig.  3).  Restoration  and 
rehabilitation  efforts,  mostly  implemented  by  federal  and  state  agencies,  attempt  to  mitigate  these 
effects,  but  often  with  mixed  results  (West,  1999;  Arkle  and  others,  2014). 
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Figure  1.  Map  showing  study  areas  (red  stars),  which  were  located  in  the  jurisdiction  of  the  Burns  Field  Office,  in 
Harney  County,  Oregon.  Base  map  sources:  Esri,  HERE,  DeLorme,  Mapmyindia,  ©  OpenStreetMap  contributors, 
and  the  GIS  user  community. 


5 


Figure  2.  Map  showing  two  transects  (red  lines)  along  eastern  facing  slopes  in  the  Stinkingwater  Mountains  (top 
panel)  and  Pueblo  Mountains  (bottom  panel),  Oregon.  We  sampled  at  three  1 -hectare  plots  (green  dots)  at  three 
elevations  along  each  transect,  for  a  total  of  18  plots  at  each  study  area.  Base  map  sources:  Esri,  HERE,  DeLorme, 
Mapmyindia,  ©  OpenStreetMap  contributors,  and  the  CIS  user  community. 
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Figure  3.  Photographs  showing  vegetation  at  the  highest  and  lowest  sampling  plots  at  the  Stinkingwater  and 
Pueblo  Mountains  study  areas.  Top  left:  Stinkingwater  high,  bottom  left:  Stinkingwater  low,  top  right:  Pueblo  high, 
and  bottom  right:  Pueblo  low. 


Study  Design  and  Sampling  Methods 

We  sampled  insects  along  an  elevation  gradient  to  quantify  associations  between  insect 
community  structure  and  climate.  The  elevation  gradient  was  established  as  uphill  transects  that 
extended  from  valley  bottom  to  ridge  tops.  Vegetation  varies  across  elevation  in  response  to  climate 
(Whittaker  andNiering,  1975).  Thus,  vegetation  and  climate  were  confounded  in  our  design.  We 
attempted  to  account  (or  control)  for  the  relative  influence  of  vegetation  on  insects  by  including 
vegetation  measurements  in  statistical  models. 

We  used  several  criteria  to  locate  and  determine  the  suitability  of  transects  for  sampling.  All 
transects  were  chosen  from  slopes  in  eastern  Oregon  where  vegetation  was  characterized  by  sagebrush 
steppe.  We  restricted  transects  to  within  1  km  of  a  road  to  provide  access.  We  intentionally  avoided 
areas  that  had  been  burned  in  the  last  25  years  or  where  vegetation  had  been  treated  (for  example, 
seeding  or  herbicide  treatments).  We  did  not  control  for  livestock  grazing  and  all  areas  occurred  within 
BLM  grazing  allotments.  These  restrictive  criteria  limited  our  ability  to  place  transects  randomly,  but 
we  attempted  to  maintain  as  much  randomization  as  possible  in  transect  and  sample  location  decisions. 
However,  given  that  we  studied  four  replicate  transects  across  two  mountain  ranges,  we  consider  the 
inference  of  this  study  to  be  approximately  sagebrush  steppe  habitats  of  eastern  Oregon. 
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Along  each  of  the  four  transects,  we  placed  nine  plots  clustered  in  groups  of  three.  A  plot  was  delineated 
as  a  1-ha  area  and  placed  within  500  m  of  the  transect  line.  The  clustering  was  a  necessity  for  sampling 
efficiency  and  allowed  for  greater  sampling  effort  at  each  elevation  along  each  transect.  Groups  of  plots 
were  placed  so  there  was  a  minimum  of  100  m  difference  in  elevation  between  them,  which  resulted  in 
variable  Euclidean  distances  between  groups  of  plots  along  transects.  Plots  within  groups  were  placed 
so  their  elevations  would  be  as  similar  as  possible,  but  the  centroids  (where  insect  samples  were  placed) 
were  always  a  minimum  of  100  m  apart. 

Insect  Sampling 

Insect  sampling  was  designed  to  quantify  the  insect  community.  We  sampled  biweekly  from 
mid-May  through  August  at  all  four  transects  in  2012  and  from  June  through  August  at  the  northern 
transect  at  the  Pueblo  Mountains  study  area  and  the  southern  transect  at  the  Stinkingwater  Mountains 
study  area  in  2013  (table  1).  We  sampled  monthly  in  September  and  October  in  2012  only. 

We  sampled  insects  using  a  protocol  described  by  Lowe  and  others  (2010).  Five  pitfall  traps 
were  placed  five  meters  from  the  center  of  each  plot  at  compass  bearings  of  36°,  108°,  180°,  252°,  and 
324°.  Pitfall  traps  were  filled  approximately  25  percent  with  low  toxicity  antifreeze.  Pitfall  traps 
contained  no  olfactory  or  visual  baits  and  therefore  were  considered  passive  traps.  Insects  were  expected 
to  fall  into  these  traps  while  walking  on  the  ground.  We  placed  one  yellow-  and  one  blue-colored 
Japanese  beetle  flight  trap  (Great  Lakes  IPM,  Inc.,  http;//www. greatlakesipm.com/)  at  each  plot,  10 
meters  from  the  center.  The  placement  of  the  first  trap  was  randomly  chosen  and  the  second  was  placed 
180°  from  the  first  (figs.  4  and  5).  The  traps  contained  no  olfactory  baits.  The  colors  of  these  traps  were 
intended  to  mimic  flowers  and  attract  flying  insects.  Blue  and  yellow  traps  are  often  used  to  mimic  the 
color  of  flowers  within  a  sampling  area,  though  there  is  evidence  that  color  of  an  insect’s  primary  host 
flower  does  not  affect  color  preference  for  traps  (Roulston  and  others,  2007;  Wilson  and  others,  2008). 
Additionally,  there  is  evidence  that  capture  rates  are  similar  between  colors  in  the  Great  Basin  Desert, 
but  that  composition  of  insect  samples  from  the  different  colors  are  different  (Wilson  and  others,  2008). 
Traps  were  left  open  for  five  consecutive  days  and  nights. 


Table  1.  Sampling  design  for  2012  and  2013,  Pueblo  and  Stinkingwater  Mountains,  Oregon. 

[Two  of  four  transects  from  2012  were  resampled  in  2013] 


Year 

Transects 

Number 
of  1-ha 
plots 

Number 
of  pitfall 
traps 

Number 
of  flight 
traps 

Number  of 
sampling 
events 

Total 

traps 

2012 

4 

36 

180 

72 

8 

2,016 

2013 

2 

18 

90 

36 

5 

630 

8 


Figure  4.  Diagram  of  pitfall  and  flight  traps  located  in  1-ha  plots.  Two  flight  traps,  one  yellow  and  one  blue,  were 
placed  1 0  meters  from  the  center  of  each  plot  and  oriented  1 80°  from  each  other.  The  outer  square  represents  the 
boundary  of  a  1  -ha  sample  plot  with  1 00  m  per  side. 


Figure  5.  Photographs  showing  pitfall  traps  (left)  with  low  toxicity  antifreeze  and  Japanese  beetle  flight  traps  (right 
top  and  bottom)  with  insecticide.  These  traps  were  used  to  capture  insects  at  each  plot,  Pueblo  and  Stinkingwater 
Mountains,  Oregon. 
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Vegetation  Sampling 

Within  each  1-ha  plot,  we  randomly  placed  six  vegetation  sampling  quadrats  with  the  constraint 
that  each  quadrat  had  to  be  at  least  20  m  from  any  other.  Each  quadrat  measured  approximately  1.5  x 
1.5  m  of  ground  area.  We  first  created  a  list  of  all  species  present  in  the  quadrat  and  then  photographed 
the  quadrat  using  a  Canon  Powershot  A590  IS  digital  camera  fixed  to  a  2  meter  monopod  (fig.  6).  This 
is  the  height  recommended  by  Booth  and  others  (2006)  for  use  of  this  technique  in  sagebrush  habitats. 

In  the  field,  we  also  recorded  the  height  of  several  functional  groups  of  vegetation;  shrubs,  native  forbs, 
native  bunch  grasses,  and  non-native  annual  grasses  that  occurred  within  each  quadrat. 

In  the  laboratory,  we  measured  percent  cover  by  species  in  each  photo-quadrat  using 
SamplePoint  software  (Booth  and  others,  2006),  which  is  a  grid-point-intercept  method  (Pilliod  and 
Arkle,  2013).  We  used  the  species  list  for  each  quadrat  that  had  been  created  in  the  field  for  reference. 
An  oversampling  effort  conducted  by  Pilliod  and  Arkle  (2013)  determined  that  six  photo-quadrats  were 
needed  to  accurately  represent  the  spatial  heterogeneity  of  vegetative  ground  cover  in  a  1-ha  sagebrush 
or  grassland  plot.  However,  earlier  work  also  had  shown  that  bunch  grasses  and  shrubs  tend  to  be 
underrepresented  in  quadrats. 

Therefore,  we  measured  the  density  and  cover  of  native  bunch  grasses  and  shrubs  using  point- 
centered  quarter  method  (hereafter  point-quarter,  fig.  6).  We  sampled  using  point-quarter  from  the 
center  of  each  of  the  six  photo-quadrats  within  each  1-ha  plot.  Within  each  of  the  four  quadrants  around 
each  sampling  point,  we  measured  the  distance  from  the  sampling  point  to  the  canopy  center  of  the 
nearest:  (1)  mature  native  bunchgrass  (canopy-line  intercept  >15  cm),  and  (2)  shrub  (canopy-line 
intercept  >10  cm),  for  a  maximum  of  4  bunchgrasses  and  4  shrubs  per  sample  point  (maximum  of  36 
bunch  grasses  and  36  shrubs  per  1-ha  plot).  Search  distances  from  sampling  points  to  the  nearest 
bunchgrass,  shrub,  or  juniper  were  limited  to  20  m  to  avoid  measuring  an  individual  plant  from  multiple 
sample  points.  We  also  recorded  the  species  and  canopy-line  intercept  of  each  plant.  Canopy-line 
intercept  was  measured  through  the  same  portion  of  the  plant  used  to  measure  the  distance  from  the 
sample  point  to  the  plant  (on  the  axis  from  the  sample  point  to  the  canopy  center).  The  intercept  value 
was  recorded  as  the  number  of  centimeters  between  the  point  on  the  meter  tape  first  intersected  by  any 
part  of  the  plant  and  the  last  point  on  the  meter  tape  intersected  by  any  part  of  the  plant  after  turning  the 
meter  tape  such  that  its  width  was  perpendicular  to  the  ground. 


Figure  6.  Photographs  showing  (left)  one  of  multiple  nadir  photographs  of  vegetation  sampling  quadrats  used  with 
SamplePoint  software  (Booth  and  others,  2006)  to  quantify  the  percent  cover  of  vegetation  and  abiotic  habitat 
characteristics  at  each  plot;  and  (right)  point-centered  quarter  method  used  to  sample  and  quantify  the  percent 
cover  of  native  bunchgrasses  and  forbs  at  each  plot. 
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Weather  and  Climate  Data 

We  measured  local  ground  temperature  by  placing  IButton®  data  loggers  (Maxim  Integrated 
Products,  Sunnyvale,  CA,  USA)  at  each  elevation  along  each  transect  (fig.  7).  We  installed  weather 
stations  at  the  mid  elevation  plots  of  each  transect  during  the  sampling  season  of  2012  and  at  high,  mid 
and  low  elevation  plots  during  the  2013  sampling  season  (fig.  7).  Weather  stations  were  installed 
approximately  one  meter  above  the  ground  and  measured  air  temperature,  solar  radiation,  and  soil 
moisture.  Weather  variables  for  plots  without  weather  stations  in  2012  were  estimated  using  iButton 
data  and  data  from  weather  stations  used  in  2013.  We  calculated  growing  degree  days  from  measured 
air  temperature  data  using  equation  1 . 


GDD  = 


Tmax  +  Tmin 
2 


10°C 


(1) 


where 


GDD  =  Growing  Degree  Days, 

Tmax  =  maximum  temperature,  and 

Tmin  =  minimum  temperature. 


Figure  7.  Photographs  showing  examples  of  how  weather  stations  were  placed  at  the  mid-elevation  plots  at  each 
transect  in  2012  and  at  every  elevation  along  each  transect  in  2013,  Stinkingwater  and  Pueblo  Mountains,  Oregon. 
Cords  and  seals  were  covered  with  duct  tape  to  provide  extra  protection  from  weather  and  wildlife.  Right 
photograph  also  shows  the  iButton®  ground  temperature  recorder  covered  by  a  solar  shield  (white  half-round  PVC). 
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Thirty-year  averages  of  climate  variables  were  collected  from  daily  DayMet  data  (Thornton  and 
others,  2014)  and  a  nearby  National  Oceanic  and  Atmospheric  Administration  weather  station  (Bald 
Mountain,  Oregon;  Latitude;  43°  33'  20"  Longitude;  118°  24'  15",  https;//www.ncdc. noaa.gov/data- 
access/land-based-station-data  https;//www. ncdc.noaa.gov/data-access/land-based-station-data)  and 
adjusted  to  estimate  local  climate  using  the  2013  iButton®  and  weather  station  data.  This  provided 
estimates  of  30-year  climate  variables  for  each  elevation  group  (three  plots)  along  each  transect.  These 
data  were  used  to  derive  biologically  meaningful  climate  variables  using  the  Bioclim  model 
(www.worldclim.org/bioclim).  For  example,  we  measured  the  30-year  average  of  the  daily  maximum 
temperature  in  the  warmest  annual  quarter  (that  is,  summer)  and  the  minimum  daily  temperature  in  the 
coldest  annual  quarter  (that  is,  winter).  The  months  included  in  these  quarters  were  variable  among 
plots,  depending  on  temperature  patterns  at  those  plots.  These  variables  provide  good  estimates  of  the 
temperature  stressors  that  insects  have  to  tolerate  within  the  sampling  plots  over  a  period  of  thirty  years. 
We  also  measured  the  30-year  average  of  cumulative  precipitation  during  the  wettest  and  driest  quarters 
at  each  plot.  These  variables  provide  good  estimates  of  long-term  seasonal  water  availability  at  the 
sampling  plots.  Future  sections  of  this  report  refer  to  weather  and  climate.  In  this  study,  we  use  the  word 
“weather”  in  reference  to  measured  (ground  and  air)  temperature  and  precipitation  variables  during  the 
sampling  period.  We  use  the  word  “climate”  in  reference  to  the  long-term  (in  this  case  30  years;  1982- 
2012)  pattern  of  air  temperature  and  precipitation  at  the  study  areas. 

Data  Management 

Percent  cover  of  vegetation  and  abiotic  environmental  characteristics  were  summarized  to  the 
plot  level  and  then  averaged  for  each  elevation.  Species  and  abiotic  characteristics  were  grouped  into 
seven  functional  groups;  soil,  rock,  bunchgrasses,  biological  crust  and  moss,  forbs  (including  non¬ 
native),  invasive  grasses,  and  shrubs.  These  data  were  derived  from  the  grid-point  intercept,  except  for 
percent  cover  of  bunchgrasses  and  shrubs,  which  were  calculated  from  point-quarter  sampling. 

Insects  were  sorted  and  identified  to  the  level  of  family.  Family-level  diversity  has  been  shown 
to  closely  parallel  species-level  diversity  in  many  taxa  including  woody  plants  (Balmford  and  others, 
1996),  ferns,  bats,  passerine  birds  (Williams  and  Gaston,  1994),  aquatic  macroinvertebrates  (Heino  and 
Soininen,  2007),  and  terrestrial  invertebrates  (Williams  and  Gaston,  1994;  Cagnolo  and  others,  2002; 
Riggins  and  others,  2009).  Restricting  identification  to  this  relatively  high  taxonomic  level  allowed  us  to 
sample  and  identify  many  more  organisms  than  if  they  had  been  identified  to  species  or  genus, 
consequently  allowing  us  to  broaden  our  temporal  sampling  and  our  area  of  inference. 
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Insects  were  analyzed  at  three  levels  of  organization:  community,  feeding  guild,  and  order 
within  feeding  guild.  At  the  finest  level  of  organization,  orders  within  feeding  guild,  we  analyzed  only 
orders  within  the  nectarivore  guild  because  nectarivorous  insects  are  important  pollinators  in  this  system 
and,  therefore,  we  had  an  interest  in  their  specific  climate  associations.  We  calculated  the  detected  insect 
family  richness  (number  of  families  captured)  and  insect  abundance  at  each  level  of  organization  at  each 
study  area  in  2012  and  2013.  We  estimated  insect  family  richness  using  the  Chao  1  richness  estimator 
(Chao  and  lost,  2012)  using  Estimates  software  (Colwell,  2013)  and  estimated  diversity  and  evenness 
using  Shannon’s  diversity  (eq.  2)  and  evenness  (eq.  3).  We  expected  that  insect  abundance  would  be 
most  strongly  influenced  by  interannual  changes  in  weather,  which  would  in  turn  affect  detectability  of 
rare  families  and  skew  richness.  To  accommodate  this  possible  sampling  bias,  we  included  Shannon’s 
diversity  and  evenness,  as  well  as  richness,  in  our  insect-weather  dataset.  Indices  of  diversity  and 
evenness  account  for  both  richness  and  abundance  of  families  within  the  community.  We  used  the  Chao 
1  richness  estimator  to  estimate  true  richness  based  on  detected  richness  and  detection  probability  (Chao 
and  lost,  2012)  to  account  for  detection  bias  associated  with  interannual  changes  in  abundance 
(detectability)  in  our  insect-climate  dataset. 

=  (2) 

i=\ 

EH=H/lnS  (3) 

where 

H'  =  Shannon’s  diversity  index, 

S  =  total  number  of  families  in  the  community. 

Pi  =  proportion  of  S  made  up  of  the  zth  species,  and 

Eh  =  evenness. 
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Data  Analysis 

In  Section  I,  Assessment  of  Sampling  Design,  we  compared  climate,  weather,  vegetation,  and 
insect  community  composition  among  study  areas  and  transects  within  study  areas.  To  make  these 
comparisons,  we  used  data  that  were  collected  at  plots  from  each  study  area  that  overlapped  in 
elevation.  We  grouped  plots  from  Stinkingwater  and  Pueblo  Mountains  study  areas  into  elevation 
classes  and  compared  plots  from  study  areas  within  elevation  classes  (fig.  8).  We  included  plots  from 
classes  2  (elevation  1,200  to  1,399  m)  and  3  (elevation  1,400  to  1,599  m),  because  these  classes  occurred 
at  both  study  areas.  We  included  vegetation  measurements  from  2012  and  2013  in  this  analysis,  to 
include  possible  interannual  variability  in  vegetation  composition.  We  used  a  general  relativization  to 
scale  variables  with  different  units  of  measurement  to  comparable  unitless  values  for  multi-response 
permutation  procedure  (MRPP).  We  used  MRPP  to  compare  samples  collected  at  different  study  areas, 
transects,  and  sampling  events.  We  used  nested  ANOVA’s  to  compare  each  weather  variable  between 
study  areas  and  among  transects  separately. 


Elevation  4 


Elevation  3 

Elevation  2 


Figure  8.  Graph  showing  sampling  plots  from  Stinkingwater  and  Pueblo  Mountains  study  areas  grouped  into  four 
elevation  classes.  The  terms  North  and  South  represent  northern  (N)  and  southern  (S)  transects  within  each  study 
area. 
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In  Section  II,  Insect  Community  Composition,  we  compare  30-year  averages  of  climate  variables 
and  interannual  weather  fluxes  from  2012  and  2013  to  insect  community  composition  along  elevation 
gradients  at  the  Stinkingwater  and  Pueblo  Mountains  study  areas.  We  used  non-parametric 
multiplicative  regression  (NPMR)  to  measure  relationships  between  30-year  averages  of  climate 
variables  and  elevation  at  both  study  areas.  For  all  of  our  NMRP  analyses,  the  neighborhood  size  was 
held  at  5  percent  of  the  sample  size.  We  included  predictors  that  improved  model  fit  to  the  response  data 
(represented  by  xR  )  by  at  least  1  percent.  We  used  PerMANOVA  to  compare  insect  community  and 
guild  richness  among  elevations.  We  explored  associations  between  insect  richness  and  long-term 
climate  variables  and  interannual  weather  fluxes  using  NPMR  and  our  Chao  1  estimates  of  richness.  For 
Section  II,  Insect  Community  Composition,  questions  1  and  2,  richness  was  calculated  using  the  entire 
2012  dataset,  including  insects  captured  in  flight  and  pitfall  traps.  For  question  3,  richness,  diversity, 
and  evenness  were  calculated  using  insects  collected  in  flight  traps  in  2012  and  2013  from  June  through 
August.  Insects  captured  in  pitfall  traps  and  insects  captured  after  August  in  2012  were  excluded  from 
analysis  in  question  3  because  sampling  was  not  conducted  after  August  in  2013  and  pitfall  samples 
were  not  collected  at  the  Pueblo  Mountains  study  area  in  2013.  The  differences  in  insect  community 
richness  between  samples  collected  in  2012  and  2013  were  compared  using  MRPP.  These  differences 
were  associated  with  the  change  in  weather  variables  between  2012  and  2013  using  NPMR.  Samples 
from  all  sampling  events  for  each  analysis  were  pooled,  and  seasonal  differences  in  richness  were  not 
considered  in  this  section. 

In  Section  III,  Insect  Phenology,  we  compare  seasonal  insect  phenology  in  2012  to  measured 
weather  variables  from  2012.  For  question  1,  we  used  blocked  MRPP  to  measure  relationships  between 
weather  and  elevation  within  sampling  events  and  to  compare  the  number  of  insect  families  that 
emerged  at  each  elevation.  We  also  used  MRPP  to  compare  the  timing  of  first  and  last  captures  of 
families  among  elevations.  We  used  NPMR  to  measure  associations  between  timing  of  first  and  last 
capture  of  insect  families  with  measured  weather.  For  question  2,  we  used  non-metric  multi-dimensional 
scaling  (NMS)  and  MRPP  to  compare  insect  community  and  guild  compositions  among  events  and 
elevations.  We  used  PerMANOVA  to  compare  abundance  of  insects  within  the  community  and  guilds 
among  events  and  elevations.  We  used  NPMR  to  measure  associations  of  insects  within  the  community 
and  guilds  with  measured  weather  variables. 
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Section  I.  Assessment  of  Sampling  Design 

We  examined  the  differences  in  climate,  weather,  vegetation,  and  insects  between  the 
Stinkingwater  and  Pueblo  Mountains  study  areas  and  among  the  four  transects  to  assess  potential  biases 
and  variation  associated  with  our  sampling  design.  We  addressed  four  primary  research  questions: 

•  Question  1:  Does  climate  vary  between  study  areas  or  between  transects  within  study  areas 
across  similar  elevations? 

•  Question  2:  Does  weather  vary  between  study  areas  or  between  transects  within  study  areas 
across  similar  elevations? 

•  Question  3:  Does  vegetation  vary  between  study  areas  or  between  transects  within  study  areas 
across  similar  elevations? 

•  Question  4:  Does  the  composition  of  insect  communities  vary  between  study  areas  or  between 
transects  within  study  areas  across  similar  elevations? 

Key  Findings 

Question  1:  Does  climate  vary  between  study  areas  or  between  transects  within  study  areas  across 
similar  elevations? 

a)  We  included  30-year  averages  of  maximum,  minimum,  and  median  daily  air  temperature  and 
precipitation  as  our  multivariate  response  in  our  MRPP  analysis  of  climate.  We  found  no 
significant  difference  in  the  30-year  averages  of  climate  between  the  study  areas  (T=-1.12, 
A=0.09,  p=0.12)  or  between  transects  within  study  areas  (Pueblo:  T=0.66,  A=-0.1 1,  p=0.73; 
Stinkingwater:  T=0.84,  A=-0.14,  p=0.79).  More  negative  values  of  T  indicate  stronger 
separation  among  groups  and  more  positive  values  of  A  represent  more  homogenous 
composition  within  groups.  Values  of  p  indicate  the  likelihood  of  a  type  II  error  and  are  used  to 
determine  statistical  significance  of  analyses. 

b)  When  we  compared  each  weather  variable  separately,  we  found  that  air  temperature  (maximum, 
minimum,  and  median)  was  significantly  higher  at  the  Pueblo  study  area  (Fi  7=16.38,  p=0.02; 
Fij=69.13,  p=0.001;  and  Fij=28.54,  p=0.006,  respectively  where  F  values  indicate  how  far  the 
data  are  dispersed  from  the  mean;  higher  F  values  indicate  larger  dispersion  and  greater 
differences  among  groups).  Precipitation  was  significantly  higher  at  the  Stinkingwater  study  area 
(Fij=39.58  p=0.003).  We  found  no  significant  differences  between  transects  within  the  plots  for 
any  of  the  measured  weather  variables. 

Question  2:  Does  weather  vary  between  study  areas  or  between  transects  within  study  areas  across 
similar  elevations  in  2012  and  2013? 

1 .  We  included  daily  average  air  temperature,  daily  average  ground  temperature,  maximum 

humidity,  and  maximum  solar  radiation  in  our  MRPP  analysis  of  weather.  When  we  compared 
study  areas  within  sampling  events  at  plots  with  overlapping  elevation,  we  found  significant 
differences  in  weather  (T=-79.88,  A=0.73,  p<0.001)  at  the  study  areas.  Pairwise  comparisons 
indicate  weather  varied  between  study  areas  through  time  (that  is,  among  sampling  events)  (table 
2). 
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Table  2.  All  pairwise  comparisons  of  weather  between  study  areas  at  different  sampling  events  were  significant  at 
p<0.001 ,  Stinkingwater  and  Pueblo  Mountain  study  areas. 


[T:  represents  the  degree  of  difference  between  groups,  in  this  case  study  areas.  A  more  negative  value  of  T  indicates 
stronger  differences  between  groups.  A:  represents  the  homogeneity  within  groups.  Larger  values  of  A  represent  more 
homogenous  groups] 


Year 

Event 

T 

A 

2012 

Early  June 

-20.28 

0.41 

Late  June 

-21.12 

0.48 

Early  July 

-22.01 

0.50 

Late  July 

-20.56 

0.43 

Early  August 

-12.79 

0.21 

Late  August 

-6.84 

0.11 

September 

-17.58 

0.33 

October 

-19.88 

0.43 

2013 

Late  May 

-10.65 

0.59 

Early  June 

-7.60 

0.33 

Late  June 

-5.83 

0.26 

Early  July 

-7.26 

0.32 

Late  July 

-7.99 

0.38 

Early  August 

-6.62 

0.30 

2.  When  we  eompared  eaeh  measured  weather  variable  separately  between  study  areas  within 
events,  we  found  signifieant  differences  in  all  of  the  variables  (average  daily  ground 
temperature:  F27, 236=6 1.56,  p<0.001;  average  daily  air  temperature:  F27, 236=340.26,  p<0.001; 
daily  maximum  humidity:  F27, 236=125. 82,  p<0.001;  and  daily  maximum  solar  radiation: 

F27, 236=208. 01,  p<0.001).  Pairwise  comparisons  of  least  squared  means  indicate  that  air 
temperature  was  significantly  higher  at  the  Pueblo  study  area  during  the  entire  season  in  2012 
and  in  the  early  season  in  2013.  Ground  temperature,  however,  was  only  significantly  different 
between  study  areas  during  one  sampling  event.  The  ground  temperature  was  significantly 
higher  at  the  Pueblo  study  area  in  early  June  2013  (Stinkingwater:  24.37°C,  Pueblo:  30.78°C, 
p=0.002).  The  Pueblo  study  area  also  had  significantly  lower  humidity  and  higher  solar  radiation 
during  most  sampling  periods.  This  pattern  was  flipped  once,  in  early  June  2013,  when  the 
humidity  at  the  Stinkingwater  study  area  dropped  below  that  at  the  Pueblo  study  area  (table  3). 
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Table  3.  Pairwise  comparisons  of  least  squares  means  at  all  sampling  events,  Stinkingwater  and  Pueblo  Mountain  study  areas,  Oregon,  2012  and 
201 3.  t=  test  statistic  that  compares  the  observed  data  to  what  would  be  expected  under  a  null  hypothesis,  p=  the  likelihood  of  a  type  II  error.  These 
are  used  to  determine  statistical  significance  of  analyses. 


Air  temperature 
(degrees  Celsius) 

Maximum  humidity 
(percent) 

Solar  radiation 
(Watts  per  square  meter) 

Year 

Event 

Pueblo 

Stinkingwater 

t 

P 

Pueblo 

Stinkingwater 

t 

P 

Pueblo 

Stinkingwater 

t 

P 

2012 

Early  June 

10.79 

10.45 

-0.66 

1.0000 

72.40 

89.68 

8.85 

<0.0001 

937.92 

784.91 

-17.05 

<0.0001 

Late  June 

21.11 

15.91 

-9.55 

<0.0001 

41.50 

65.03 

11.87 

<0.0001 

953.28 

893.17 

-6.69 

<0.0001 

Early  July 

25.01 

19.26 

-10.22 

<0.0001 

39.42 

60.65 

10.89 

<0.0001 

943.96 

885.90 

-6.47 

<0.0001 

Late  July 

25.29 

21.32 

-6.72 

<0.0001 

38.28 

59.22 

10.81 

<0.0001 

882.60 

827.04 

-6.19 

<0.0001 

Early  August 

25.97 

22.96 

-4.95 

0.0004 

33.35 

45.49 

5.10 

0.0002 

927.29 

846.30 

-8.92 

<0.0001 

Late  August 

26.33 

24.6 

2.36 

1.0000 

43.11 

45.32 

0.93 

1.0000 

814.48 

719.11 

-9.32 

<0.0001 

September 

19.85 

16.88 

5.12 

0.0002 

30.14 

42.36 

7.55 

<0.0001 

814.15 

719.11 

-10.59 

<0.0001 

October 

11.13 

8.25 

4.81 

0.0008 

64.02 

79.95 

8.15 

<0.0001 

572.10 

452.57 

-13.32 

<0.0001 

2013 

Late  May 

12.98 

9.03 

5.31 

<0.0001 

59.98 

78.00 

7.15 

<0.0001 

1,064.37 

840.37 

-17.65 

<0.0001 

Early  June 

24.07 

20.08 

5.30 

<0.0001 

41.25 

76.82 

1.66 

1.0000 

934.77 

903.29 

-2.48 

1.0000 

Late  June 

16.39 

17.61 

-1.64 

1.0000 

63.41 

77.22 

4.10 

0.0189 

896.11 

872.72 

-5.39 

<0.0001 

Early  July 

26.59 

26.05 

0.72 

1.0000 

57.54 

44.61 

-0.75 

0.1000 

913.33 

915.83 

0.20 

1.0000 

Late  July 

28.14 

26.84 

1.74 

1.0000 

23.29 

34.37 

3.29 

0.4101 

974.02 

874.78 

-7.82 

<0.0001 

Early  August 

23.15 

21.28 

2.51 

1.0000 

32.86 

44.26 

3.39 

0.2934 

903.49 

835.53 

-5.35 

0.0070 
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c)  Within  study  areas,  we  compared  the  same  weather  variables  between  sampling  transects 
and  events  from  2012.  We  did  not  include  2013  in  this  analysis  because  we  sampled  only 
one  transect  per  study  area  in  2013.  We  found  significant  differences  at  the  Stinkingwater 
and  Pueblo  study  areas  (T=-43.49,  A=0.96,  p<0.0001;  and  T=-38.84,  A=0.71,  p<0.001, 
respectively).  Pairwise  comparisons  indicate  that  transects  within  study  areas 
experienced  different  weather  for  every  sampling  event  (table  4). 

d)  We  found  no  significant  differences  in  air  temperature  and  few  significant  differences  in 
ground  temperature  between  transects  at  the  Pueblo  study  area.  Ground  temperatures 
were  significantly  higher  at  the  northern  transect  in  early  July  (least  squares  means  at 
northern  transect:  26.45°C  and  southern  transect:  18.10°C,  p=0.0002)  and  significantly 
higher  at  the  southern  transect  in  October  (least  squares  means  at  northern  transect: 
20.84°C  and  southern  transect:  26.32°C,  p<0.0001).  Similarly,  there  were  few  differences 
in  air  or  ground  temperature  between  transects  at  the  Stinkingwater  study  area.  In  late 
June  and  early  August,  air  temperature  was  higher  at  the  northern  transect  (least  squares 
means  at  northern  transect:  17.1  TC  and  southern  transect:  14.72°C  p=0.04;  and  at 
northern  transect:  24.35°C  southern  transect:  21.57°C,  p=0.004,  respectively).  In  early 
June  ground  temperatures  were  higher  at  the  southern  transect  (least  squares  means  at 
northern  transect:  15.64°C  and  southern  transect:  17.5 TC,  p=0.02).  However,  we  found 
significant  differences  in  solar  radiation  between  transects  in  many  sampling  events  at 
both  study  areas  (table  5).  We  found  significant  differences  in  solar  radiation  between 
transects  at  both  study  areas  during  every  sampling  event,  though  these  differences  were 
not  consistent  between  transects  within  the  study  areas. 


Table  4.  Multi-response  permutation  procedure  comparisons  of  weather  between  transects  at  every  sampling 
event  at  the  Pueblo  and  Stinkingwater  Mountain  study  areas,  Oregon,  2012. 


Stinkingwater  Mountains 

Pueblo  Mountains 

Event 

T 

A 

P 

T 

A 

P 

Early  June 

-11.6 

0.91 

<0.001 

-5.19 

0.24 

<0.001 

Late  June 

-10.29 

0.57 

<0.001 

-6.28 

0.25 

<0.001 

Early  July 

-11.5 

0.83 

0.007 

-5.03 

0.19 

<0.001 

Late  July 

-10.5 

0.59 

0.006 

-5.41 

0.21 

<0.001 

Early  August 

-11.36 

0.78 

<0.001 

-4.43 

0.18 

0.002 

Late  August 

-8.26 

0.39 

0.100 

-6.52 

0.26 

<0.001 

September 

-11.34 

0.75 

0.012 

-5.06 

0.22 

<0.001 

October 

-11.61 

0.94 

<0.001 

-4.88 

0.22 

<0.001 
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Table  5.  Solar  radiation  was  significantly  different  between  transects  within  Stinkingwater  and  Pueblo  Mountain 
study  areas,  Oregon,  2012.  Transects  at  the  Stinkingwater  Mountains  study  area  were  19  km  apart  and  transects 
at  the  Pueblo  Mountains  study  area  were  9  km  apart. 


Solar  Radiation 
[watts  per  square  meter) 

Study  Area 

Event 

Northern  transect 

Southern  transect 

t 

P 

Early  June 

818.13 

757.69 

166.27 

<0.0001 

Late  June 

897.52 

888.82 

21.77 

<0.0001 

Early  July 

859.8 

911.99 

-130.62 

<0.0001 

Stinkingwater 

Late  July 

834.78 

819.3 

38.74 

<0.0001 

Mountains 

Early  August 

806.03 

886.56 

-206.19 

<0.0001 

Late  August 

736.1 

725.65 

26.04 

<0.0001 

September 

694.15 

744.1 

-124.91 

<0.0001 

October 

399.4 

505.78 

-266.32 

<0.0001 

Early  June 

935.43 

1064.37 

12.43 

<0.0001 

Late  June 

997.4 

934.77 

-220.58 

<0.0001 

Early  July 

963.12 

896.11 

-95.89 

<0.0001 

Pueblo 

Late  July 

914.78 

913.33 

-161.08 

<0.0001 

Mountains 

Early  August 

942.72 

974.02 

-77.2 

<0.0001 

Late  August 

860 

903.49 

-227.85 

<0.0001 

September 

853.32 

940.4 

-196.03 

<0.0001 

October 

595.83 

909.15 

-118.91 

<0.0001 

Question  3:  Does  vegetation  vary  between  study  areas  or  between  transeets  within  study  areas  aeross 
similar  elevations? 

a)  We  found  significant  differences  in  percent  cover  of  vegetation  functional  groups 

between  study  areas  within  elevation  class  2  (T=-9.42,  A=0.16,  p<0.001),  but  not  class  3 
(T=-0.49,  A=0.02,  p=0.30).  This  indicates  that  vegetation  composition  was  different  at 
the  study  areas  at  lower  elevations,  but  was  more  similar  at  higher  elevations.  We  found 
no  significant  differences  in  vegetation  composition  between  transects  within  study 
areas. 
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Question  4:  Does  the  eomposition  of  insect  communities  vary  between  study  areas  or  between  transects 
within  study  areas  across  similar  elevations? 

a)  When  we  compared  insect  communities  between  study  areas  using  PerMANOVA,  we 
found  that  the  insect  community  composition  at  the  Pueblo  and  Stinkingwater  study 
areas  were  significantly  different  (Fi_  192=2. 69,  p=0.05).  Event  also  was  significant  in  this 
analysis  (F14, 192=8.02,  p<0.001).  This  analysis  does  not  account  for  differences  in 
elevation  of  the  sampling  plots  at  the  study  areas.  Therefore,  the  analysis  was  run  a 
second  time  with  a  subset  of  samples  from  each  study  area  that  overlapped  in  elevation 
(fig.  8).  This  second  analysis  used  fewer  samples,  but  accounted  for  the  repeated 
measures  (events)  and  controlled  for  elevation.  We  found  a  nearly  significant  difference 
between  study  areas  for  this  analysis  (Fi_  i28=2.26,  p=0.06). 


b)  We  found  that  transects  were  significantly  different  from  one  another  (Fa^  i96=3.0, 
p<0.001;  fig.  9).  Pairwise  comparisons  indicated  that  this  significant  relationship  was 
driven  by  the  southern  transect  at  the  Pueblo  study  area,  which  was  significantly  different 
from  the  north  transect  at  the  same  study  area  (t=l  .73,  p=0.006)  and  the  southern  transect 
at  the  Stinkingwater  study  area  (t=2.05,  p<0.001).  The  significant  relationship  between 
transects  was  maintained  when  the  analysis  was  run  with  a  subset  of  plots  from  similar 
elevations  (F3ji2=2.51,  p=0.007). 
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Figure  9.  Non-metric  multidimensional  scaling  analysis  shows  that,  regardless  of  elevation,  insect  community 
composition  at  plots  from  the  Pueblo  Mountains  study  area  were  more  similar  to  each  other  than  to  that  at  the  plots 
from  the  Stinkingwater  Mountains  (Stink)  study  area.  Insect  communities  at  plots  within  a  transect  (indicated  by 
PuebloN,  PuebloS,  StinkN,  StinkS)  were  also  more  similar  to  each  other  than  to  other  transects,  with  a  few 
exceptions.  Plots  from  the  highest  elevation  cluster  at  the  Pueblo  study  area  were  most  similar  to  plots  from  the 
lowest  elevation  cluster  at  the  Stinkingwater  Mountains  study  area. 
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Interpretation 

Although  climate  was  not  significantly  different  between  the  study  areas,  the  Pueblo  study  area 
is  consistently  hotter  and  drier  than  the  Stinkingwater  study  area,  and  this  trend  persisted  during  2012 
and  2013.  Vegetation  was  significantly  different  between  the  study  areas  at  the  low  elevation  plots. 
These  differences  may  be  influenced  by  differences  in  land  use  at  the  low  plots  because  low  elevation 
areas  are  often  more  accessible  and,  therefore,  more  heavily  used  for  cattle  grazing  and  recreation.  It 
also  could  be  affected  by  slope,  which  has  been  shown  to  change  heatload  and  direct  radiation 
(McCune,  2007).  Insect  community  composition  differed  both  between  study  areas  and  between 
transects  within  study  areas.  These  differences  represent  variability  in  the  system  associated  with  biotic 
and  abiotic  habitat  characters. 
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Section  II.  Insect  Community  Composition 

We  examined  relationships  between  insects  and  climate  at  three  levels  of  biological 
organization:  entire  insect  community,  insects  organized  by  guilds,  and  insect  pollinators  grouped  by 
order.  We  addressed  three  primary  research  questions: 

•  Question  1:  How  do  climate  and  vegetation  vary  along  elevation  gradients? 

•  Question  2:  How  do  insect  communities  vary  in  relation  to  climate  and  vegetation? 

•  Question  3:  How  do  interannual  fluctuations  in  weather  affect  insect  composition  and  how  do 

these  relationships  compare  to  long-term  relationships  between  insects  and  climate? 

We  predicted  that  plots  at  high  elevations  would  have  cooler  temperatures  and  high  annual 
precipitation.  This  would  directly  cause  differences  in  the  composition  of  vegetation  and  insect 
communities  at  the  plots.  In  addition  to  climate,  we  predicted  that  vegetation  composition  would 
directly  affect  the  composition  of  insect  communities  (fig.  10).  We  expected  that  long-term  climate 
variables  would  have  the  greatest  influence  on  the  richness  of  insect  communities  and  that  short-term 
weather  variables  would  more  strongly  affect  diversity  and  evenness  indices,  because  they  are 
dependent  on  abundance  as  well  as  richness. 

We  captured  59,517  insects  from  176  families  and  10  orders  at  the  Pueblo  study  area  and 
1 12,305  insects  from  185  families  and  1 1  orders  at  the  Stinkingwater  study  area  in  2012  and  2013.  Of 
all  the  individuals  captured  at  the  Stinkingwater  study  area,  77,688  were  from  the  family  Cecidomyiidae 
(Diptera,  gall  gnats). 
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Figure  10.  A  depiction  of  the  predicted  relationships  among  elevation,  climate,  weather,  vegetation,  and  insects. 
We  predicted  that  elevation  would  affect  both  climate  and  weather.  We  expected  that  vegetation  and  insect 
community  compositions  would  be  associated  with  changes  in  climate  and  weather  and  that  insect  communities 
would  also  vary  in  association  with  changes  in  vegetation  composition. 
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Key  Findings 

Question  1:  How  do  climate  and  vegetation  vary  along  elevation  gradients? 

The  30-year  average  air  temperature  during  the  aetive  season  (May-September)  deereased  with 
elevation  and  the  30-year  average  eumulative  annual  preeipitation  inereased  with  elevation  at  our  study 
plots  (fig.  1 1).  Differenees  in  air  temperature  are  most  likely  to  affeet  inseets  during  the  aetive  season, 
beeause  they  spend  their  inaetive  season  in  diapause,  whieh  allows  them  to  avoid  extreme  temperatures. 
However,  eumulative  annual  preeipitation  is  a  good  measure  of  response  to  moisture,  beeause  inaetive 
season  preeipitation  direetly  affeets  aetive  season  moisture  availability  through  snowpaek.  Most  elimate 
variables  measured  were  assoeiated  with  elevation  in  the  NPMR  analysis,  although  the  strengths  of  the 
relationships  varied  (table  6).  The  strongest  relationships  were  between  the  standard  deviation 
(hereinafter,  variability)  in  average  daily  air  temperatures  (daily  maximum,  daily  minimum,  and  median 
air  temperatures)  and  elevation,  with  xR  values  greater  than  0.90.  The  fit  of  a  model  to  the  dataset  ean 
be  deseribed  by  xR^  and  ean  be  interpreted  similarly  to  R^  values  in  a  linear  regression.  However,  the 
xR  is  ealeulated  using  an  iterative  proeess  of  removing  one  datapoint  at  a  time  and  estimating  the  value 
of  that  point  based  on  the  other  data  in  the  set.  Beeause  of  the  differenee  in  how  xR^  is  ealeulated,  very 
weak  models  may  have  negative  xR  values,  whieh  ean  be  interpreted  as  0  or  no  fit  (MeCune,  2011). 
The  variability  in  average  daily  air  temperature  was  lower  for  plots  at  higher  elevations.  The 
relationship  between  the  variability  of  preeipitation  and  elevation  was  stronger  than  that  between 
average  preeipitation  and  elevation.  The  variability  of  preeipitation  reeeived  among  years  was  greater 
for  plots  at  higher  elevations  than  for  relatively  lower  plots. 
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Figure  11.  Linear  regression  analysis  showing  that  plots  at  high  elevations  had  low  30-year  average  air 
temperatures  and  high  30-year  average  precipitation,  Stinkingwater  and  Pueblo  Mountain  study  areas,  Oregon. 
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Table  6.  Associations  between  climate  and  elevation  indicate  that  the  variability  of  air  temperature  among  years, 
rather  than  maximum,  median,  or  minimum  air  temperatures,  has  the  strongest  association  with  elevation  at  the 
Stinkingwater  and  Pueblo  Mountain  study  areas,  Oregon. 

[High  xR^  values  indicate  a  strong  non-linear  relationship  between  elevation  and  the  climate  variable] 


Climate  Variable 

xR^ 

Variance  of  maximum  daily  air  temperature 

0.98 

Variance  of  median  daily  air  temperature 

0.97 

Variance  of  minimum  daily  air  temperature 

0.91 

Maximum  air  temperature  warmest  quarter 

0.68 

Variance  of  minimum  air  temperature  coldest  quarter 

0.56 

Variance  of  precipitation  wettest  quarter 

0.48 

Variance  of  precipitation  driest  quarter 

0.45 

Maximum  daily  air  temperature 

0.45 

Variance  of  maximum  air  temp  warmest  quarter 

0.25 

Variance  of  annual  precipitation 

0.23 

Median  daily  air  temperature 

0.20 

Annual  precipitation 

0.03 

Minimum  air  temperature  coldest  quarter 

0.01 

Precipitation  wettest  quarter 

<0.01 

a)  Percent  cover  of  forbs,  shrubs  and  invasive  grasses  were  not  associated  with  elevation.  However, 
bunchgrasses  had  higher  percent  cover  at  higher  elevations  (xR  =0.28). 

b)  Although  most  vegetation  variables  were  not  associated  with  elevation,  they  were  associated 
with  climate,  particularly  the  variability  of  climate  variables.  Percent  cover  of  shrubs  was  higher 
at  plots  with  more  variable  precipitation  among  years  (xR  =0.4621).  This  suggests  that  shrubs 
need  at  least  some  years  with  higher  amounts  of  precipitation,  though  they  are  more  likely  to 
persist  through  drought  years  than  other  vegetation  taxa  with  shallower  rooting  systems,  such  as 
bunchgrasses.  Percent  cover  of  bunchgrasses  was  higher  at  plots  with  less  variable  high  air 
temperatures  in  the  summer  and  less  variable  minimum  air  temperatures  in  the  winter 
(xR^=0.5482).  Invasive  annual  grasses  had  the  highest  percent  cover  at  plots  that  had  moderate 
minimum  air  temperatures,  and  low  variability  in  precipitation  during  the  driest  season,  which 
tended  to  be  the  drier  sites  at  lower  elevations  (xR^=0.7024).  Finally,  forbs  had  the  highest 
percent  cover  at  plots  that  had  the  highest  variability  in  precipitation  in  the  wettest  and  driest 
quarters  (xR^=0.4621).  Plots  with  more  variable  climate  patterns  during  the  growing  season 
contained  more  native  vegetation  and  plots  that  were  consistently  drier  contained  more  invasive 
annual  grasses. 


25 


Question  2:  How  do  insect  communities  vary  in  relation  to  climate  and  vegetation? 

a)  We  found  no  significant  difference  in  insect  community  composition  among  elevation 
groups  at  the  Pueblo  study  area  (F=  1.10,  p=0.320)  but  there  was  a  significant  difference 
at  the  Stinkingwater  study  area  (F=  3.41,  p<0.001).  Pairwise  comparisons  indicate  that 
the  significant  relationship  was  driven  by  differences  between  the  highest  set  of  plots 
(elevation  1,859  m  to  1,880  m)  with  the  lowest  two  sets  of  plots  (1,282  m  to  1,307  m). 

b)  Insect  community  richness  at  the  sampling  plots  was  not  associated  with  any  of  the 
climate  variables  measured  for  this  study.  However,  guild-level  analyses  of  insect 
richness  indicated  that  the  variability  in  climate  variables  among  years  was  often  more 
strongly  related  to  richness  of  the  guild  than  the  values  of  the  measured  variables 
themselves  (table  7).  Four  feeding  guilds  were  associated  with  climate  variables: 
detritivores,  nectarivores,  parasites,  and  predators  (table  7).  Herbivores  were  the  only 
guild  that  was  not  associated  with  any  measured  environmental  variable.  Predator 
richness  was  higher  in  plots  where  the  minimum  winter  air  temperature  was  relatively 
moderate.  Detritivore  richness  was  highest  at  plots  with  relatively  moderate  variability  in 
summer  high  temperatures.  Nectarivore  richness  was  higher  at  plots  with  more  variability 
in  minimum  winter  air  temperature  and  more  variability  in  precipitation.  Finally, 
parasites  had  the  highest  richness  at  plots  with  relatively  moderate  variability  in  cold 
season  temperatures  and  where  precipitation  was  highest  during  the  driest  (and  warmest) 
season. 

c)  The  richness  of  detritivores  and  predators  were  the  only  guild  measurements  associated 
with  vegetation  (table  7).  Detritivores  were  positively  associated  with  shrubs  and 
predators  were  negatively  associated  with  non-native  annual  grasses. 


Table  7.  Guild-level  associations  between  richness  of  insect  families  and  environmental  variables  indicate  that 
insects  are  associated  with  climate  and  vegetation  at  the  Stinkingwater  and  Pueblo  Mountain  study  areas,  Oregon. 

[High  xR^  values  indicate  a  strong  non-linear  relationship  between  feeding  guild  and  the  climate  variable] 


Guild 

xR2 

Environmental  Predictor  1 

Environmental  Predictor  2 

Detritivore 

0.0984 

Variance  of  annual  maximum  air 
temperature 

Percent  cover  shrubs 

Nectarivore 

0.2249 

Variance  of  daily  minimum  air 
temperature 

Variance  of  annual  precipitation 

Parasite 

0.3222 

Variance  of  median  air  temperature  in  the 
coldest  quarter 

Precipitation  in  the  driest  quarter 

Predator 

0.2125 

Percent  cover  annual  grass 

Minimum  air  temperature  in  the  coldest 
quarter 
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d)  We  captured  pollinator  families  (nectarivores)  from  three  orders:  Diptera  (flies), 

Hymenoptera  (bees),  and  Lepidoptera  (butterflies  and  moths).  The  richness  of  flies  was 
associated  with  percent  cover  of  forbs  and  the  variability  of  annual  extreme  cold  air 
temperatures  (xR  =0.761 1).  Plots  with  higher  percent  cover  of  forbs  had  greater 
variability  in  the  number  of  nectarivorous  fly  families  captured.  Plots  with  high  forb 
cover  contained  two  nectarivorous  fly  families,  bombylidae  (bee  flies)  and  syrphidae 
(hover  flies).  Plots  with  less  variability  in  air  temperatures  in  the  winter  had  higher 
richness  than  plots  with  more  moderate  variability  in  winter  air  temperatures  due  to  the 
addition  of  one  family,  Conopidae  (thick-headed  flies).  Plots  with  higher  variability  in 
their  minimum  winter  air  temperatures  and  less  annual  grass  supported  higher  richness  of 
bees  (xR  =0.5843).  The  plots  with  higher  annual  grass  cover  also  were  less  variable  in 
the  richness  of  bees  than  plots  with  lower  annual  grass  cover.  The  richness  of  butterflies 
and  moths  was  slightly  lower  and  less  variable  at  plots  with  higher  percent  cover  of 
forbs,  possibly  resulting  from  increased  competition  at  these  plots.  Richness  of  this  group 
was  also  lower  at  plots  with  more  variability  in  minimum  cold  air  temperatures  among 
years  (xR^=0.1 127). 

Question  3:  How  do  interannual  fluctuations  in  weather  affect  insect  composition  and  how  do  these 
relationships  compare  to  long-term  relationships  between  insects  and  climate? 

Flying  Insects 

a)  Variation  in  insect  composition  between  2012  and  2013  was  measured  at  all  study  areas 
for  flying  insects  (through  flight  traps),  but  only  Stinkingwater  for  crawling  insects 
(through  pitfall  traps).  We  found  that  flying  insect  community  composition  was 
significantly  different  between  sampling  years  at  both  the  Stinkingwater  (T=-8.56, 
A=0.30,  p<0.001)  and  Pueblo  (T=-8.51,  A=0.29,  p<0.001)  study  areas.  The  insect 
communities  at  both  study  areas  changed  similarly,  with  a  reduction  in  richness, 
diversity,  and  evenness  from  2012  to  2013.  MRPP  analysis  indicated  that,  although 
insect  community  composition  at  the  study  areas  were  significantly  different  (T=-2.93, 
A=0.05,  p=0.02),  the  magnitude  and  direction  of  the  change  in  the  community 
compositions  between  years  at  the  study  areas  was  not  significantly  different  (T=-0.82, 
A=0.03,  p=0.18).  This  indicates  that  flying  insect  community  composition  shifted  in 
response  to  environmental  factors  similarly  at  both  study  areas.  We  found  no  significant 
difference  in  the  composition  of  crawling  insects  between  years  at  the  Stinkingwater 
study  area  (T=-0.49,  A=0.02,  p=0.23). 

b)  We  found  associations  between  interannual  changes  in  weather  and  richness,  diversity, 
and  evenness  of  the  flying  insect  community  (table  8).  Richness  had  the  lowest  xR  value 
(0.2880),  indicating  that  relationships  between  changes  in  richness  and  weather  were 
more  difficult  to  predict  than  those  between  changes  in  diversity  or  evenness  and  weather 
at  this  level  of  organization. 
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Table  8.  Community-level  analyses  of  changes  in  insect  richness,  diversity,  and  evenness  with  changes  in 
environmental  variables  between  2012  and  2013  indicate  that  the  community  is  affected  by  interannual  variability  in 
weather  and  vegetation  at  Stinkingwater  and  Pueblo  Mountain  study  areas,  Oregon. 

[High  xR^  values  indicate  a  strong  non-linear  relationship  between  the  insect  community  statistic  and  the  climate  variable] 


Statistic 

xR2 

Environmental  Predictor  1 

Environmental  Predictor  2 

Richness 

0.288 

Study  area 

Median  daily  ground  temperature 

Diversity 

0.5546 

Maximum  air  temperature  in  the 
coldest  quarter 

Percent  cover  annual  grass 

Evenness 

0.4689 

Minimum  air  temperature  in  the  coldest 
quarter 

Percent  cover  annual  grass 

c)  Decreased  richness,  diversity,  and  evenness  of  the  flying  insect  community  from  2012  to 
2013  were  associated  with  changes  in  vegetation  and  weather  including  increased  non¬ 
native  annual  grass  cover,  decreased  median  ground  temperatures  in  the  active  season, 
and  variable  maximum  air  temperatures  in  the  winter  (decreased  at  some  plots  and 
increased  at  some  plots).  The  variability  in  diversity  and  evenness  of  insects  among  plots 
that  had  stable  or  increasing  maximum  winter  air  temperatures  between  2012  and  2013 
was  much  smaller  than  the  variability  of  plots  that  had  cooler  cold  season  air 
temperatures  in  2013.  This  suggests  that  milder  (i.e.  warmer)  winters  may  benefit  some 
families,  but  not  others.  Finally,  plots  that  had  similar  ground  temperatures  in  the  active 
season  for  both  years,  or  that  were  warmer  in  2013  than  2012,  retained  more  richness 
than  plots  that  had  cooler  ground  temperatures  in  2013. 

d)  Changes  in  richness  within  guilds  was  more  strongly  associated  with  changes  in  weather 
than  with  climate  (tables  7  and  9),  except  for  nectarivores.  Changes  in  richness  of 
nectarivores  were  relatively  poorly  associated  with  changes  in  the  measured  weather 
variables.  Evenness  and  diversity  of  predators,  nectarivores,  and  herbivores  were 
strongly  associated  with  changes  in  weather,  but  this  relationship  was  not  as  strong  for 
detritivores  or  parasites  (table  9).  Differences  in  the  strength  of  relationships  between 
insect  guilds  and  the  weather  and  climate  variables  with  which  they  were  associated 
indicate  that  guild-level  analyses  reveal  trends  that  are  specific  to  life  history  traits 
associated  with  feeding  guilds. 

e)  All  guilds  were  associated  with  some  measure  of  warm  season  temperatures,  whether  air 
or  ground  temperature  (table  9).  Three  guilds,  herbivore,  nectarivore,  and  parasite,  were 
also  sensitive  to  temperature  in  the  cold  season,  when  insects  are  in  diapause.  For  most 
guilds,  plots  with  the  greatest  change  in  summer  air  temperatures  experienced  the 
greatest  decreases  in  richness,  diversity,  and  evenness  whereas  plots  with  the  greatest 
change  in  ground  temperature  (during  the  active  season)  experienced  the  smallest  losses 
in  these  measures.  Most  plots  had  colder  winter  air  temperatures  in  2013  than  in  2012, 
but  the  plots  that  cooled  less  or  (the  few)  that  were  warmer  in  2013  retained  higher 
richness  and  diversity  of  herbivores,  nectarivores  and  parasites.  Contrary  to  what  we 
expected,  richness  associated  more  strongly  with  changes  in  environmental  variables  (for 
example,  vegetation  and  weather)  between  2012  and  2013  than  richness  associated  with 
climate  variables. 
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Pollinators 


f)  We  found  that  nectarivores  were  associated  with  weather  variables.  The  associations 
between  flies  and  weather  variables  were  not  as  strong  as  those  between  flies  and 
climate.  The  temperature  variables  with  which  flies  associated  were  warm  season 
variables,  as  opposed  to  the  cold  season  variables  with  which  they  were  associated  in  the 
climate  analysis  (table  10).  However,  their  richness  was  higher  at  warmer  plots  in  the 
warm  season  and  lower  at  the  colder  plots  in  the  cold  season,  consistently  indicating  a 
sensitivity  to  cold  air  temperatures.  Bees,  moths,  and  butterflies  associated  with  similar 
climate  and  weather  variables.  As  in  the  climate  analysis,  the  association  between  the 
richness  of  moths  and  butterflies  and  weather  variables  was  not  very  strong,  although 
their  diversity  and  evenness  were  more  strongly  associated  with  weather. 

g)  Bees,  butterflies,  and  moths  responded  to  interannual  changes  in  vegetation.  The 
evenness  of  butterflies  and  moths  and  the  richness  of  bees  were  retained  at  higher  levels 
at  plots  that  lost  fewer  bunchgrasses  between  2012  and  2013,  although  the  evenness  of 
butterflies  and  moths  was  variable  among  plots  where  there  was  no  loss  of  bunchgrasses. 
The  diversity  of  bees  was  associated  with  the  percent  cover  of  shrubs.  Most  plots  gained 
bee  diversity  between  years,  and  the  plots  that  gained  the  most  diversity  also  gained  the 
most  shrub  cover.  Butterflies  and  moths  were  the  only  order  within  the  nectarivore  guild 
to  associate  with  forbs.  Although  moths  and  butterflies  decreased  overall  between  years, 
plots  with  greater  increases  in  forb  cover  between  years  had  smaller  decreases  in  richness 
of  moths  and  butterflies. 

h)  Flies  stood  out  from  the  other  two  orders  because  they  associated  only  with  changes  in 
summer  air  temperature  variables,  whereas  bees,  butterflies,  and  moths  also  associated 
with  changes  in  winter  air  temperatures  and  vegetation  variables.  Overall,  most  plots  lost 
richness,  diversity,  and  evenness  of  flies  between  2012  and  2013,  but  plots  with  warmer 
summer  air  temperatures  (minimum,  mid,  and  maximum)  retained  more  richness  and 
diversity  than  plots  with  cooler  summer  air  temperatures.  Bees,  butterflies,  and  moths 
were  sensitive  to  air  temperatures  during  the  winter.  Diversity  and  evenness  of  bees 
increased  between  sampling  years  in  response  to  warmer  winters.  Detected  richness  of 
bees,  however,  had  U-shaped  distribution  where  plots  with  little  or  no  change  in  winter 
air  temperatures  lost  the  most  richness  and  plots  with  more  extreme  cold  or  warm  winters 
retained  more  richness.  These  results  are  counter-intuitive,  but  it  is  possible  that  lower 
winter  temperatures  indicate  more  consistent  snow  insulation  for  individuals  buried  right 
below  the  surface  of  the  soil.  Warmer  winters  may  simply  be  more  survivable  for  over¬ 
wintering  pupae  even  without  snow  insulation.  Richness,  diversity,  and  evenness  of 
butterflies  and  moths  all  decreased  at  most  plots  between  sampling  years  in  response  to 
lower  median  winter  temperatures  in  2013.  However,  plots  that  cooled  less  between 
years  retained  higher  detected  richness,  diversity  and  evenness  of  moths  and  butterflies. 


29 


Table  9.  Guild-level  analyses  of  the  relationship  between  interannual  changes  in  insect  richness,  diversity,  and  evenness  and  interannual  changes 
in  environmental  variables  at  the  Stinkingwater  and  Pueblo  Mountain  study  areas,  Oregon. 

[High  xR^  values  indicate  a  strong  non-linear  relationship  between  the  guild  and  the  climate  variable] 


Guild 

Statistic 

xR2 

Environmental  Predictor  1 

Environmental  Predictor  2 

Detritivore 

Richness 

0.5992 

Annual  daily  max  ground  temperature 

Daily  median  ground  temperature 

Diversity 

0.3943 

Elevation 

Daily  minimum  air  temperature  in  the  warmest  quarter 

Evenness 

0.3538 

Elevation 

Percent  cover  bunchgrass 

Herbivore 

Richness 

0.8076 

Study  area 

Daily  minimum  air  temperature  in  the  coldest  quarter 

Diversity 

0.7929 

Daily  maximum  ground  temperature 

Percent  cover  shrubs 

Evenness 

0.6979 

Annual  grass 

Percent  cover  shrubs 

Nectarivore 

Richness 

0.2299 

Daily  max  air  temperature  in  the  warmest  quarter 

Percent  cover  forbs 

Diversity 

0.7701 

Elevation 

Daily  median  air  temperature  in  the  coldest  quarter 

Evenness 

0.797 

Elevation 

Daily  minimum  air  temperature  in  the  warmest  quarter 

Parasite 

Richness 

0.7319 

Daily  median  air  temperature  in  the  coldest  quarter 

Percent  cover  annual  grass 

Diversity 

0.1404 

Daily  minimum  air  temperature  in  the  warmest  quarter 

Percent  cover  bunchgrass 

Evenness 

0.0691 

Daily  median  air  temperature  in  the  coldest  quarter 

Daily  median  air  temperature  in  the  warmest  quarter 

Predator 

Richness 

0.6776 

Daily  median  air  temperature  in  the  warmest  quarter 

Daily  median  ground  temperature 

Diversity 

0.6228 

Daily  max  ground  temperature 

Percent  cover  bunchgrass 

Evenness 

0.7759 

Daily  minimum  air  temperature  in  the  warmest  quarter 

Percent  cover  bunchgrass 
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Table  10.  Analyses  of  orders  within  the  nectarivore  guild  indicate  that  pollinators  from  different  orders  are  sensitive  to  changes  in  different  weather 
and  vegetation  variables.  Stinkingwater  and  Pueblo  Mountain  study  areas,  Oregon. 

[High  xR^  values  indicate  a  strong  non-linear  relationship  between  the  order  and  the  climate  variable] 


Order 

Statistic 

xR2 

Environmental  Predictor  1 

Environmental  Predictor  2 

Diptera 

Richness 

0.106 

Daily  maximum  air  temperature  in  the  warmest  quarter 

Daily  median  air  temperature  in  the  warmest  quarter 

Diversity 

0.281 

Daily  minimum  air  temperature  in  the  warmest  quarter 

Daily  maximum  ground  temperature 

Evenness 

0.2908 

Study  area 

Daily  median  ground  temperature 

Hymenoptera 

Richness 

1 

Daily  maximum  air  temperature  in  the  coldest  quarter 

Percent  cover  bunchgrass 

Diversity 

0.5212 

Daily  median  air  temperature  in  the  coldest  quarter 

Percent  cover  shrub 

Evenness 

0.7609 

Daily  median  air  temperature  in  the  coldest  quarter 

Daily  minimum  ground  temperature 

Lepidoptera 

Richness 

0.1966 

Daily  median  air  temperature  in  the  coldest  quarter 

Percent  cover  forbs 

Diversity 

0.5112 

Daily  median  air  temperature  in  the  coldest  quarter 

Daily  minimum  ground  temperature 

Evenness 

0.5594 

Daily  median  air  temperature  in  the  coldest  quarter 

Percent  cover  bunchgrass 
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Interpretation 

Variability  in  Climate,  Weather,  and  Habitat  Affected  Richness,  Diversity,  and  Evenness  of  Insect 
Communities 

We  found  evidence  that  climate  affected  insect  community  composition  and  structure  at  the 
scale  of  this  study.  We  found  that  air  temperature  variability,  rather  than  long-term  average  air 
temperature,  was  most  often  associated  with  insect  community  composition.  The  relationship  between 
precipitation  and  elevation  was  weak,  indicating  that  there  may  not  have  been  enough  variability  in 
elevation  to  detect  relationships  between  precipitation  and  insects.  We  only  detected  a  relationship 
between  parasites  and  the  30-year  average  of  precipitation.  It  is  possible  that  a  greater  change  in 
elevation  between  sampling  plots  would  have  revealed  more  relationships  because  insect  populations 
may  contain  enough  genetic  plasticity  to  accommodate  elevational  differences  in  temperature  through 
altered  body  color  and  size,  wing  formation,  fecundity,  or  timing  of  life  history  traits  (Hodkinson, 
2005).  Changes  in  body  size  along  elevation  gradients  may  decrease  metabolic  costs  (Buckley  and 
others,  2013). 

At  the  scale  of  this  study,  insect  communities  appeared  to  be  more  strongly  influenced  by 
weather  than  climate.  Our  results  suggest  that  interannual  variation  in  weather,  especially  air 
temperature,  limit  the  distribution  of  insects.  Insect  species  are  often  limited  by  the  amount  of  thermal 
energy  they  need  to  complete  development  within  their  active  season  (Bird  and  Hodkinson,  1999; 
Hodkinson,  2005).  Unpredictable  swings  in  air  temperature  may  limit  the  distribution  of  some  insect 
groups  by  limiting  the  ability  of  juveniles  to  successfully  survive  to  adulthood  in  years  with  low  air 
temperatures. 

Effects  on  Pollinators 

Insect  pollinator  assemblages  have  been  shown  to  change  along  elevation  gradients  in  response 
to  changes  in  climate  and  vegetation.  Species  compositions  tend  to  move  from  bee-dominated  at  lower 
elevations  and  fly-dominated  at  higher  elevations  (Kearns,  1992;  Arnold  and  others,  2009),  at  least 
partially  in  response  to  differences  in  humidity  (Devoto  and  others,  2005). 

Although  bees  have  been  shown  to  dominate  relatively  dry  habitat  compared  to  flies,  our 
analyses  indicate  that  bees  have  lower  family  richness  where  invasive  annual  grasses  are  more 
abundant,  which  tend  to  be  lower  elevation,  dry  plots  (Knutson  and  others,  2014).  Fly  pollinators  at 
higher  elevations  may  compete  with  bees,  which  could  limit  the  range  of  some  bees  along  elevation 
gradients  (Kearns,  1992).  Increased  intensity  and  frequency  of  disturbance  events,  such  as  fire,  have 
been  shown  to  increase  the  spread  of  invasive  annual  grasses  at  drier,  lower  elevation  areas  of  the  Great 
Basin  (Balch  and  others,  2013).  This  disturbance  may  reduce  and  fragment  habitat  and  reduce  habitat 
quality  for  bee  pollinators  in  this  system. 

Invasive  Annual  Grasses  Affected  the  Distribution  of  Insects  Among  Plots 

Land  use  change  and  invasive  species  interact  with  changing  climate  to  increase  the  effects  of 
natural  stressors  on  communities  and  change  disturbance  regimes.  These  changes  may  be  as  important 
as  or  more  important  than  the  direct  effect  of  change  in  climate  variables  in  predicting  ecosystem 
response  to  climate  change  (Dunne  and  others,  2004).  We  found  evidence  of  the  effect  of  invasive 
annual  grasses  on  insects  at  every  level  of  analysis.  By  replacing  native  vegetation,  non-native  annual 
grasses  may  exacerbate  habitat  fragmentation  that  could  lead  to  reduced  foraging  success  for  bees  or 
hunting  success  for  predators. 
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Section  III.  Insect  Phenology 

The  goal  of  this  part  of  the  study  was  to  document  phenological  shifts  in  insect  communities  in 
response  to  temperature  along  elevation  gradients.  We  examined  the  following  questions: 

•  Question  1:  How  is  the  timing  of  insect  emergence  or  diapause  related  to  variability  in  weather 
along  an  elevation  gradient  at  the  beginning  and  end  of  the  active  season? 

•  Question  2:  How  does  insect  abundance  vary  through  time  in  relation  to  environmental 
characteristics,  such  as  air  and  ground  temperature,  growing  degree  days,  relative  humidity,  and 
solar  radiation? 

We  predicted  that  insects  emerge  later  at  higher  elevations  due  to  lower  springtime  temperatures 
and,  therefore,  have  shorter  annual  activity  periods  than  insects  at  lower  elevations. 

Key  Findings 

Question  1:  How  is  the  timing  of  insect  emergence  or  diapause  related  to  variability  in  weather  along 
an  elevation  gradient  at  the  beginning  and  end  of  the  active  season? 

a)  Weather  was  significantly  variable  among  elevations  within  sampling  events  as  indicated  by 
blocked  MRPP  analyses  (T=-13.041,  A=0.021,  p<0.0001).  We  measured  cumulative  growing 
degree  days  (formula  1)  and  maximum  humidity  at  the  plots  and  found  that  they  were 
significantly  different  among  elevations  within  sampling  events  (T=-9.491  A=0.016  p<0.0001 
and  T=-2 1.401  A=0.199  p<0.0001,  respectively). 

b)  We  compared  the  number  of  families  that  were  captured  for  the  first  time  at  each  elevation.  We 
found  significantly  more  families  were  captured  for  the  first  time  at  lower  elevations  in  both 
years  (2012:  T=-40.36,  A=0.08,  p<0.001;  2013:  T=-9.91,  A=0.04,  p<0.001;  fig.l2).  Similarly, 
significantly  more  families  were  captured  for  the  last  time  at  lower  elevations  in  2012  (T=-6.37, 
A=0.01,  p<0.001;  2013,  fig.  12).  However,  we  found  that,  throughout  the  entire  sampling  season 
in  2012,  more  insect  families  were  captured  at  higher  elevations  at  the  Pueblo  (F5, 17=20.37, 
p<0.001)  and  Stinkingwater  (F5j7=4.17,  p=0.02)  study  areas. 
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Figure  12.  Bubble  graph  showing  the  number  of  insect  families  captured  for  the  first  time  or  the  last  time  at  each 
sampling  elevation  and  at  what  time  during  the  season  they  were  captured,  Stinkingwater  and  Pueblo  Mountains, 
Oregon.  The  size  of  the  bubbles  represents  the  number  of  families  captured  and  are  scaled  relative  to  each  other. 
As  a  point  of  reference,  the  bubble  in  the  top  left  corner  of  each  panel  is  labeled  with  the  number  of  families  that  it 
represents.  Most  families  were  captured  for  the  first  time  early  in  the  sampling  season  at  relatively  low  elevations. 
Though  the  relationship  was  not  as  dramatic,  most  families  were  captured  for  the  last  time  at  relatively  low 
elevations  late  in  the  season,  often  in  September  or  October,  2012. 
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c)  We  found  that  timing  of  first  capture  was  significantly  different  among  elevations  for  every 
feeding  guild  in  2012  (table  11).  Time  of  first  capture  was  later  at  higher  elevations  for  every 
guild  except  detritivores,  which  emerged  earlier  at  higher  elevations  in  2012  (fig.  13).  Timing  of 
last  capture  was  significantly  different  among  elevations  for  herbivores,  predators  and  parasites 
in  2012  (table  1 1).  We  did  not  measure  last  capture  in  2013,  because  our  sampling  was  truncated 
to  focus  on  the  early  part  of  the  season. 

d)  We  found  strong  correlations  between  the  timing  of  first  capture  of  insects  from  feeding  guilds 
with  weather  variables  in  2012  (table  12;  fig.  13).  Minimum  ground  temperature,  humidity,  solar 
radiation,  and  growing  degree  days  were  the  most  important  variables  that  we  measured  in  2012. 


Table  11.  Multi-response  permutation  procedure  analysis  indicates  differences  in  timing  of  first  and  last  captures  of 
families  within  feeding  guilds  among  elevations  in  2012. 

[T :  represents  the  separation  between  groups,  the  more  negative  T,  the  greater  the  separation.  A:  represents  the  homogeneity 
within  groups,  a  larger  value  of  A  represents  a  group  that  is  more  homogenous,  p=  the  likelihood  of  a  type  II  error] 


First  Capture 

2012 

Guild 

T 

A 

P 

n 

Herbivore 

-15.91 

0.096 

<0.0001 

461 

Detritivore 

-4.03 

0.II4 

0.0007 

lOI 

Nectarivore 

-9.77 

0.092 

<0.0001 

295 

Predator 

-3.84 

0 

0.0012 

237 

Parasite 

-5.85 

O.I 

<0.0001 

261 

Last  Capture 

2012 

Guild 

T 

A 

P 

n 

Herbivore 

-3.102 

0.036 

0.0089 

192 

Detritivore 

1.306 

0.08 

0.133 

33 

Nectarivore 

0.2349 

0.006 

0.348 

96 

Predator 

1. 173 

0.064 

0.029 

83 

Parasite 

-1.3309 

0.039 

0.034 

II8 
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Figure  13.  Average  timing  of  first  capture  for  all  families  grouped  by  feeding  guilds  was  significantly  different 
among  elevations  in  2012  (large  windows,  table  11),  Stinkingwater  and  Pueblo  Mountains,  Oregon.  Weather 
variables  measured  along  elevation  gradients  indicated  that  timing  of  first  capture  was  strongly  associated  with 
measured  weather  variables  for  all  feeding  guilds  (small  windows,  table  12).  The  Y-axis  of  the  small  window  plots 
(showing  weather  associations)  represents  the  number  of  families  within  the  guild  that  were  captured  for  the  first 
time  in  any  sample. 
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Figure  13 —Continued 


Table  12.  Timing  of  first  and  last  capture  of  insect  families  within  guilds  was  strongly  associated  with  measured  weather  variables  at  the 
Stinkingwater  and  Pueblo  Mountain  study  areas,  Oregon.  In  this  analysis,  some  guilds  were  best  explained  by  three  predictor  variables  while  others 
were  best  explained  by  two  predictor  variables. 

[High  xR^  values  indicate  a  strong  non-linear  relationship  between  feeding  guild  and  the  climate  variable] 


First  Capture 

2012 

Guild 

xR^ 

Environmental  predictor  1 

Environmental  predictor  2 

Environmental  predictor  3 

Detritivore 

0.6509 

Daily  minimum  ground  temperature 

Daily  maximum  humidity 

Growing  degree  days 

Herbivore 

0.7845 

Daily  minimum  ground  temperature 

Daily  maximum  solar  radiation 

Growing  degree  days 

Nectarivore 

0.7418 

Minimum  ground  temperature  recorded 

Daily  maximum  humidity 

Growing  degree  days 

Parasite 

0.6781 

Daily  minimum  ground  temperature 

Maximum  solar  radiation  recorded 

Growing  degree  days 

Predator 

0.7551 

Event 

Maximum  solar  radiation  recorded 

Last  Capture 

2012 

Guild 

xR^ 

Environmental  predictor  1 

Environmental  predictor  2 

Environmental  predictor  3 

Detritivore 

0.0966 

Minimum  air  temperature  on  coldest  day 

Maximum  humidity  on  most  humid  day 

Herbivore 

0.6517 

Event 

Daily  maximum  ground  temperature 

Nectarivore 

0.5969 

Event 

Daily  maximum  humidity 

Parasite 

0.1961 

Daily  maximum  solar  radiation 

Minimum  day  length 

Predator 

0.6017 

Event 

Growing  degree  days 
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Question  2:  How  does  insect  abundance  vary  through  time  in  relation  to  environmental  characteristics, 
such  as  air  and  ground  temperature,  relative  humidity,  and  solar  radiation? 

a)  Measured  weather  (air  temperature,  ground  temperature,  humidity,  and  solar  radiation)  differed 
among  sampling  events,  and  among  elevations  within  events  (table  13,  p<0.001  for  all  analyses). 
In  both  years,  MRPP  values  of  T  were  more  negative  and  values  of  A  were  larger  for  samples 
grouped  by  event,  indicating  that  these  groups  were  more  separated  from  each  other  and  more 
homogenous  in  composition  than  groups  defined  by  elevation  within  events.  The  mean 
difference  between  maximum  air  temperatures  among  the  lowest  and  highest  elevations  within 
events  on  each  transect  was  8.72  ±  0.49°C  in  2012  and  6.95  ±  2.69°C  in  2013.  The  mean 
difference  between  the  highest  and  lowest  maximum  air  temperatures  at  a  single  elevation 
among  the  warmest  and  coolest  sampling  events  was  25.17  ±  1.70“C  in  2012  and  24.52  ±  2.16°C 
in  2013. 

b)  Insect  composition  varied  significantly  among  elevations  within  events  (T=-21 .364,  A=0.290, 
P>0.001)  and  among  events  (T=-49.684,  A=0.152,  p>0.001;  figs.  14  and  15)  in  2012  on  the  basis 
of  MRPP  and  NMS. 

c)  We  compared  the  timing  of  peak  abundance  for  all  families  that  were  captured  at  every  sampling 
event,  and  found  that  the  timing  of  peak  abundance  for  the  most  abundant  insect  families  was 
staggered  throughout  the  season  (fig.  15). 

d)  As  in  Section  II.  Insect  Community  Composition,  we  grouped  insect  families  into  feeding  guilds 
for  analysis.  Herbivores  and  nectarivores  were  the  most  abundant  insect  guilds  captured 
throughout  the  sampling  season  in  both  years.  Both  of  these  guilds  peaked  in  abundance  in  the 
early  part  of  the  season,  though  herbivores  also  had  a  second  larger  peak  at  the  end  of  the  season 

(fig.  16). 


Table  13.  Multi-response  permutation  procedure  analysis  was  used  to  compare  weather  among  sampling  events 
and  among  elevations  nested  within  sampling  events  at  our  plots  located  at  the  Stinkingwater  and  Pueblo 
Mountains,  Oregon.  In  2012  and  2013,  there  were  significant  differences  among  plots  measured  repeatedly 
(events)  and  plots  measured  during  the  same  events  along  an  elevation  gradient  (p<0.0001  for  all  analyses). 

[T :  represents  the  separation  between  groups,  the  more  negative  T,  the  greater  the  separation.  A:  represents  the  homogeneity 
within  groups,  a  larger  value  of  A  represents  a  group  that  is  more  homogenous,  p=  the  likelihood  of  a  type  II  error] 


Year 

T 

A 

Group 

2013 

-28.93 

0.410 

Event 

-12.81 

0.302 

Elevation  within  event 

2012 

-86.65 

0.579 

Event 

-32.06 

0.485 

Elevation  within  event 

39 


Axis  1  (42.2%) 


Figure  14.  Non-metric  multidimensional  scaling  showing  that  insect  samples  from  different  events  during  the 
sampling  season  changed  somewhat  systematically,  with  the  mean  value  for  each  month  shifting  through 
ordination  space. 
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•  Coccinellidae  3; 

Empididae  910 

• 

Platygastridae  25 
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Figure  15.  Peak  abundance  of  families  captured  during  every  sampling  event,  Stinkingwater  and  Pueblo 
Mountains,  Oregon.  The  sizes  of  the  bubbles  represent  the  relative  abundance  of  the  families  at  peak  abundance. 
The  bubbles  are  also  labeled  with  the  number  of  individuals  they  represent.  The  bubbles  are  colored  coded  by 
order. 


2012  Feeding  Guilds  through  Time 

12000 


-•—Herbivore  Detritivore  Nectarivore-*- Predator  -•-Parasite 


Figure  16.  Abundances  of  insects  by  feeding  guilds  throughout  the  sampling  season  in  2012.  All  feeding  guilds 
peaked  in  abundance  in  June  and  July  in  2012,  but  herbivores  and  nectarivores  were  by  far  the  most  frequently 
captured  insects.  Herbivores  had  a  second  peak  in  October. 
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d)  We  found  that  the  eomposition  of  insect  guilds  varied  significantly  by  sampling  event 
and  by  elevations  within  sampling  events,  indicating  that  seasonal  variability  in  weather 
at  plots  and  differences  in  weather  among  plots  associated  with  elevation  affected  the 
composition  of  insect  guilds  (table  14).  Two  guilds,  herbivores  and  predators,  varied  in 
abundance  among  elevations  within  events  (fig.  17). 

f)  We  used  NPMR  to  associate  changes  in  weather  resulting  from  event  and  elevation  to  the 
abundance  of  insect  guilds.  We  included  event  and  elevation  as  predictors  in  the  model  to 
represent  variability  that  was  not  explained  by  the  measured  weather  variables.  Minimum 
air  temperature  and  elevation  explained  32.4  percent  of  the  variance  in  the  abundance  of 
detritivores  among  samples,  the  highest  xR  value  achieved  for  any  guild  in  this  analysis 
(table  15).  Samples  from  high  elevations  with  higher  minimum  air  temperatures 
(probably  from  warmer  sampling  events)  contained  fewer  detritivores.  The  abundance  of 
predators  was  also  relatively  well  predicted  by  elevation  and  maximum  humidity  with 

3 1 .4  percent  of  the  variance  explained.  In  this  case,  plots  and  events  with  higher 
humidity  had  higher  abundances  and  plots  at  higher  elevation  had  slightly  lower 
abundances.  Parasites  were  most  strongly  associated  with  temperature,  ground  and  air 
temperatures  explained  18.4  percent  of  the  variability  in  this  group.  We  captured  slightly 
fewer  parasites  at  plots  and  events  with  higher  temperatures.  The  variability  in 
abundance  of  herbivores  and  nectarivores  was  not  well  explained  by  our  model,  6.2  and 

8.4  percent,  respectively. 

g)  Within  the  nectarivore  guild,  bees  were  the  most  abundant  pollinators,  especially  in  June 
and  early  July  (fig.  18).  In  August,  butterflies  and  moths  became  more  abundant  than 
bees,  but  their  overall  abundance  throughout  the  season  was  lower.  When  we  compared 
the  abundance  of  insect  orders  within  the  nectarivore  guild  using  PerMANOVA,  we 
found  that,  as  in  the  guild-level  analysis,  all  orders  within  the  nectarivore  guild  varied 
significantly  in  composition  among  events  and  elevations  within  events  (table  16). 

Unlike  the  guild-level  analysis,  we  found  that  the  abundance  of  all  of  the  orders  also 
changed  significantly  among  events  and  elevation  within  events  (table  16). 

Table  14.  Composition  of  insect  guilds  was  significantly  different  among  sampling  events  and  elevations  within 
events. 

[F=the  magnitude  of  differences  among  groups,  p=  the  likelihood  of  a  type  II  error] 


Elevation 
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Guild 
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P 

F 
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8.679 
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0.0002 
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0.0002 
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0.0002 
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0.0002 
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0.0002 

6.079 

0.0002 
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Figure  17.  Abundance  of  herbivores  and  predators  varied  among  elevations  within  sampling  events,  Stinkingwater 
and  Pueblo  Mountains,  Oregon.  The  elevation  with  the  highest  abundance  in  each  sampling  event  is  drawn  in  bold 
and  has  the  number  of  individuals  captured  near  or  inside  of  it. 


2012  Pollinators  by  Event 
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Figure  18.  Abundance  of  pollinators  from  early  June  to  October  2012  at  the  Stinkingwater  and  Pueblo  Mountain 
study  areas,  Oregon.  Bees  were  the  most  abundant  pollinator  in  the  early  season,  butterflies  and  moths  were  more 
commonly  collected  later  in  the  season. 
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Table  15.  Guilds  were  associated  with  weather  variables  to  varying  degrees  at  the  Stinkingwater  and  Pueblo 
Mountain  study  areas,  Oregon.  In  this  analysis,  some  guilds  were  best  explained  by  three  predictor  variables  while 
others  were  best  explained  by  two  predictor  variables.  The  abundances  of  detritivores  and  predators  were  most 
strongly  associated  with  weather,  parasites  were  moderately  associated,  and  herbivores  and  nectarivores  were 
only  weakly  associated  with  the  measured  weather  variables. 

[High  xR^  values  indicate  a  strong  non-linear  relationship  between  the  feeding  guild  and  the  climate  variable.] 


Guild 

xR2 

Variable  1 

Variable  2 

Variable  3 

Detritivore 

0.3244 

Elevation 

Minimum  air 
temperature 

Herbivore 

0.0623 

Growing  degree  days 

Maximum  solar 
radiation 

Nectarivore 

0.0841 

Elevation 

Maximum  humidity 

Parasite 

0.1842 

Variance  of  ground 
temperature 

Maximum  air 
temperature 

Minimum  air 
temperature 

Predator 

0.3135 

Elevation 

Maximum  humidity 

Table  16.  Differences  in  composition  and  abundances  in  nectarivorous  orders  of  insects  among  events  and 
elevations  within  events. 

[Orders  include  nectarivorous  members  of  bees  (Hymenoptera),  flies  (Diptera),  and  butterflies  and  moths  (Lepidoptera); 
F=the  magnitude  of  differences  among  groups,  p=  the  likelihood  of  a  type  II  error] 


Metric 

Elevation 

Event 

Composition 

Order 

F 

P 

F 

P 

Hymenoptera 

3.2 

0.0002 

16.221 

<0.0001 

Diptera 

3.352 

0.0002 

14.654 

0.0002 

Lepidoptera 

2.668 

0.0002 

15.041 

0.0002 

Abundance 

Hymenoptera 

13.077 

0.0002 

3.0428 

0.0002 

Diptera 

20.127 

0.0002 

2.7434 

0.0002 

Lepidoptera 

21.854 

<0.0001 

2.0631 

0.0002 
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h)  When  we  associated  the  abundance  of  nectarivorous  insect  orders  with  weather  variables, 
we  found  that  flies  were  the  most  strongly  associated  group,  with  21.9  percent  of  the 
variance  in  abundance  explained  by  the  variability  of  air  temperatures  and  solar  radiation 
(table  17).  Bees,  butterflies,  and  moths  were  poorly  associated  with  the  measured 
weather  variables,  with  only  3.5  and  1.8  percent  of  the  variation  in  abundance  explained. 
Although  associations  for  two  groups  were  weak,  all  three  were  associated  with  solar 
radiation.  Plots  and  events  with  higher  solar  radiation  had  fewer  pollinators. 


Table  17.  Associations  between  nectarivorous  insect  orders  and  weather  variables.  Flies  were  the  only  pollinator 
group  that  associated  strongly  with  the  measured  weather  variables  at  our  study  areas. 

[Orders  include  nectarivorous  members  of  bees  (Hymenoptera),  flies  (Diptera),  and  butterflies  and  moths  (Lepidoptera); 
High  xR^  values  indicate  a  strong  non-linear  relationship  between  elevation  and  the  climate  variable.] 


xR2 

Variable  1 

Tolerance 

Variable  2 

Tolerance 

Diptera 

0.2186 

STDAir 

3.41281 

MaxSolar 

89.75416 

Hymenoptera 

0.0349 

GDD 

1,309.807 

Max  Solar 

89.75416 

Lepidoptera 

0.0175 

STDGmd 

1.99265 

Max  Solar 

448.7708 

Interpretation 

Weather,  Elevation,  and  the  Timing  of  Emergence  and  Diapause 

We  found  evidence  that  the  emergence  and  diapause  timing  of  insects  is  dependent  on 
environmental  cues  associated  with  seasonal  and  elevational  differences  in  weather.  We  captured  more 
families  for  the  first  and  last  time  at  low  elevation  plots,  but  more  families  overall  at  high  elevation 
plots.  This  could  indicate  that  insects  from  the  same  groups  have  longer  active  seasons  at  lower 
elevations,  or  that  some  insect  taxa  may  have  migrated  seasonally  along  an  elevation  gradient  to  track 
favorable  weather  conditions.  Seasonal  migrations  of  insects  in  temperate  zones  often  involve 
movements  to  and  from  diapause  sites  (Dingle,  1982).  Ashton  and  others  (2009)  found  that  larvae  of  a 
rare  mountain  butterfly  species  modify  their  behavior  along  an  elevation  gradient  to  regulate 
temperature.  If  climate  change  patterns  continue  as  predicted  and  high  elevation  sites  continue  to  warm, 
it  is  possible  that  insect  taxa  that  regulate  temperature  through  seasonal  migration  up  and  down 
elevation  gradients  may  become  stranded  in  unsuitable  habitat  at  the  tops  of  warming  elevation  dines, 
especially  during  the  hottest  and  driest  times  in  their  active  season. 
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Variability  in  Guild  Abundance  Throughout  the  Active  Season 

We  found  that  the  most  abundant  families  had  peak  abundanees  that  were  staggered  through  the 
aetive  season.  Seasonal  ehanges  in  temperature  and  preeipitation  were  greater  than  ehanges  among 
elevations  within  the  same  sampling  period.  These  staggered  peaks  in  abundanee  may  be  examples  of 
adaptive  seasonality,  whieh  allows  inseet  taxa  to  avoid  temperatures  outside  of  their  thermal  toleranee 
or,  in  the  ease  of  thermally  tolerant  taxa,  to  avoid  eompetition  from  more  thermally  limited  taxa  (Powell 
and  Logan,  2005).  Univoltine  speeies,  sueh  as  most  inseets  in  temperate  regions,  are  espeeially  sensitive 
to  eontinuously  inereasing  temperatures  that  are  assoeiated  with  elimate  ehange.  Relatively  infrequent, 
extreme  (high)  temperature  events  ean  result  in  large  deereases  in  reproduetion  due  to  redueed  survival 
of  immature  individuals  (Kingsolver  and  others,  2011). 

Pollinators,  Seasonality,  and  Climate  Change 

Pollinators  are  a  group  of  partieular  interest  due  to  their  eeonomie  and  eeologieal  importanee. 

We  found  that  the  peak  abundanees  of  bees  and  butterflies,  the  most  abundant  pollinators,  were 
staggered  through  time.  Bees  from  the  family  Halietidae  were  espeeially  abundant  in  the  early  summer 
and  moths  from  the  family  Tortrieidae  were  espeeially  abundant  in  the  late  summer.  Colonial 
thermoregulation  is  eommon  in  soeial  bees,  sueh  as  European  honey  bees  (Kronenberg  and  Heller, 

1982),  but  Halietid  bees,  as  well  as  most  native  bees  in  this  region,  are  solitary.  In  this  ease,  bees  may 
avoid  foraging  (and  pollinating)  at  high  temperatures  (Heinrieh,  1975).  High  temperatures  have  also 
been  shown  to  adversely  affeet  reproduetion  in  Tortrieid  moths  (Carroll  and  Quiring,  1992).  Large 
females  generally  produeed  more  eggs  under  eontrol  temperatures,  but  had  redueed  feeundity  and 
longevity  under  eonditions  of  inereased  temperature. 

Management  Implications  and  Future  Directions 

Measurements  of  ehanges  in  elimate  and  weather  along  elevation  gradients  provide  a  unique 
opportunity  to  study  the  effeets  of  elimate  on  eommunities  while  eontrolling  for  faetors  sueh  as  genetie 
adaptation  in  or  among  populations.  These  faetors  are  eonfounded  with  elimate  effeets  in  studies  that 
measure  responses  over  large  areas  eneompassing  many  latitudes  (Hill  and  others,  2011;  Simberloff, 
2009). 

Information  derived  from  our  analyses  ean  be  applied  to  make  predietions  regarding  the  eurrent 
and  future  distribution  of  inseets  in  a  highly  disturbed  and  eeonomieally  important  semi-arid  eeosystem. 
The  analysis  teehnique  and  information  produeed  by  this  projeet  eould  be  used  in  future  studies  to  help 
prediet  habitat  quality  and  loeations  for  organisms  of  interest,  sueh  as  pollinators. 

The  plan  for  the  seeond  phase  of  this  projeet  is  to  assess  ehanges  in  speeies  ranges  and 
phenology  10  years  or  more  after  this  study  at  the  same  established  sampling  plots.  Ideally  this  would 
oeeur  in  the  years  2022-2025,  beeause  many  elimate  foreeasts  prediet  ehanges  for  the  years  2025  and 
2050.  We  eould  address  questions  sueh  as:  (1)  Have  speeies  ranges  shifted  over  time?  (2)  How  have 
eommunity  eomposition  and  speeies  abundanee  and  phenology  ehanged  over  time?  (3)  Have  ehanges 
been  uniform  throughout  elevation  gradients?  (4)  Have  some  speeies  been  affeeted  more  strongly  at 
speeifie  elevations?  This  seeond  phase  of  the  projeet  eould  provide  important  data  that  ean  be  used  to 
develop  better  speeies-elimate  models  and  predict  how  insect  communities  will  be  affected  by  continued 
climate  change. 
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